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ABSTRACT 

Ambipolar diffusion is important in redistributing magnetic flux and in damping Alfven waves in 
molecular clouds. The importance of ambipolar diffusion on a length scale £ is governed by the 
ambipolar diffusion Reynolds number, i?AD = ^/-^ad, where Iad is the characteristic length scale for 
ambipolar diffusion. The logarithmic mean of the AD Reynolds number in a sample of 15 molecular 
clumps with measured magnetic fields ()Crutcheiill999D is 17, comparable to the theoretically expected 
value. We identify several regimes of ambipolar diffusion in a turbulent medium, dependi ng on the 
ratio of the flow time to collision times between ions and neutrals; the clumps observed bv iCrutcheii 
(jl999() are all in the standard regime of ambipolar diffusion, in which the neutrals and ions are coupled 
over a flow time. We have carried out two-fluid simulations of ambipolar diffusion in isothermal, 
turbulent boxes for a range of values of -Rad- The mean Mach numbers were fixed at A^ = 3 and 
Ma = 0.67; self-gravity was not included. We study the properties of overdensities-i.e., clumps- 
in the simulation and show that the slope of the higher-mass p ortion of the c lump mass spectrum 
increases as Rad decreases, which is qualitatively consistent with iPadoan et all {2007) 's finding that 
the mass spectrum in hydrodynamic turbulence is significantly steeper than in ideal MHD turbulence. 
For a value of Rad similar to the observed value, we find a slope that is consistent with that of the 
high-mass end of the Initial Mass Function for stars. However, the value we find for the spectral 
index in our ideal MHD simulation differs from theirs, presumably because our si mulations have 
different initial conditions. This suggests that the mass spectrum of the clumps in the IPadoan et al.l 
((20 07') turbulent fragmentation model for the IMF depends on the environment, which would conflict 
with evidence for a universal IMF. In addition, we give a general discussion of how the results of 
simulations of magnetized, turbulent, isothermal boxes can be scaled to physical systems. Each 
physical process that is introduced into the simulation, such as ambipolar diffusion, introduces a 
dimensionless parameter, such as Rad, which must be fixed for the simulation, thereby reducing the 
number of scaling parameters by one. We show that the importance of self-gravity is fixed in any 
simulation of ambipolar diffusion; it is not possible to carry out a simulation in which self-gravity and 
ambipolar diffusion are varied independently unless the ionization is a free parameter. We show that 
our simulations apply to small regions in molecular clouds, generally with io < 0.4 pc and M < 25 Mq. 
A general discussion of the scaling relations for magnetized, isothermal, turbulent boxes, including 
self-gravitating systems, is given in the Appendix. 

Subject headings: Magnetic fields — MHD — ISM: magnetic fields — ISM: kinematics and dynamics — 
stars:formation 



1. INTRODUCTION 

Giant molecular clouds, threaded by magnetic fields, 
are the birth places for new stars. Since the earliest 
studies of star formation, it has been recognized that 
the magnetic flux in stars is many orders of magnitude 
less than that in th e interstellar material f rom which 
the stars originated. iMestel fc: Spitzed ([1956') suggested 
that ambipolar diffusion (AD) could resolve this prob- 
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lem by allowing magnetic flux to be redistributed during 
collapse due to the differential motion between the ion- 
ized and neutral gas. With effective shielding of high 
energy cosmic rays and radiation, the ionization fraction 
of gas inside high-d e nsity cloud cores ca n be < 10^^ (e.g. 
Caselli et al.l 119981: iBergin et al] I1999D . which renders 
AD efficient. Star formation theory based on the AD- 
regulated, qu a si-static collapse of niolecul a r clouds (e.g. 



Spitzer" '1968'; 'Nakano & Tadcmaru" '1972"; 



Mouschovia^ 
1976, 1977, 1 979; Nakano fc Nakamura 1978; Shu 198S 
Li^ano fc Shillfl98a iFiedJer fc Mouschmoasl 119921 fl993 
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naturally accounts for the enormous loss of magnetic flux 
during star formation. 

Ho we ver, both observatio ns (iZuckerman fc Evans! 
119741 iZuckerman fc Palmed 119741 ) and theory 
(jArons fc Maxl[l975l ) have long indicated that supersonic 
turbulent motions are important in molecular clouds, 
and this turbulence h a s a ma j or effect on star formation 
(iMac Low fc Klesseiil 120041: iBallesteros-Paredes et alj 
l2007HMcKee fc Ostrikeil2007[) . The kinetic energy of the 
supersonic motions is observed to be comparable to the 
magnetic energy of the clouds, so that i nolecular clouds 
are in approximate equipartition (e.g. [C rutchcr' J1999I : 
IHeiles fc TroTandl [200a iTroland fc Cnitcher. .20081) . It 
should be borne in mind that the amplitude of turbulent 
fluct uations decrease s with decr easing scale; for exam- 
ple. [Goodman et "aD ([1998) and iBarranco fc GoodmanI 
([1998') find that the NH3 lines within ~ 0.1 pc of the 
centers of the cores that they examined do not obey the 
line width-size relations as seen on the large scale (e.g . 
Larson" 198 11; 'Solo mon et al.lll987l;lHever fc Brunt' l200l. 
Mouschovias (19871 and iMvers fc Lazaria n (1998) have 
argued that AD would d amp turbulen t motions on 
small scales, and [Goodma n et al.l ()1998[ ) suggest that 
this damping could be enhanced by the low ionization 
fraction in the dense inner regions of the cores. Better 
data will enable determination of the role of ambipolar 
diffusion on the small scales relevant to the formation 
of individual stars. Existing data clearly show that 
turbulence is important on larger scales, but observa- 
tional tests of the theoretical p rediction that turbulence 
can accelerat e the rate of AD (jFatuzzo fc Adarni 120021 : 
lZweibel|[200^ will be challenging. 

Numerical simulation is an important tool in under- 
standing supersonic turbulence in magnetized MCs, but 
it is very challenging to carry out three-dimensional (3D) 
simulations that include ambipolar diffusion. The small 
ionization fraction in molecular clouds means that the 
ion inertia can be neglected. This permits a single-fluid 
treatment of ambipolar diffusion, which gives the in- 
duction equation the form of a diffusion eq uation (e.g. 
IMac Low et an[l995l : Imiffin fc Pudritzj 120081) . However, 
in this case the the stability condition for explicit codes 
requires the time step to scale as the square of the grid- 
size (Aa;^ — [Mac Low et al."1995 D, which is prohibitive at 
high resolution (e.g. iNa kamura fc Lill200"8[ ). Li, McKee, 
fc Klein (2006; hereafter LMK) developed the Heavy-Ion 
Approximation, which takes advantage of the negligible 
ion inertia in regions of very low ionization and can accel- 
erate simulations of ambipolar diffusion by large factors. 
In the Heavy- Ion Approximation, the mass- weighted ion- 
ization is increased by a factor 7?, ~ 10^ and the ion- 
neutral coupling coefficient is decreased by the same fac- 
tor, so that the momentum transfer between ions and 
neutrals is unaffecte d. Using a semi-implicit two-fluid 
scheme pro posed bv IMac Low fc SmithI ()1997[ ) (see also 
lTdthlll995l ). LMK tested the Heavy- Ion Approximation 
with several classical problems involving ambipolar dif- 
fusion and found speed-ups of order a factor 100. 

In the first astrophysical application of the Heavy-Ion 
Approximation, Li et al (2008, hereafter LMKF) stud- 
ied the statistical properties of supersonically turbulent 
systems with ambipolar diffusion. The properties of the 
turbulence were found to vary smoothly from the hydro- 
dynamic case to the ideal MHD case as the importance 



of ambipolar diffusion decreased. They found that the 
power spectra for the neutral gas properties of a strongly 
magnetized medium with strong ambipolar diffusion are 
similar to those for a weakly magnetized medium; in par- 
ticular, the power spectrum for the neutral velocity is 
close to that for Burger's turbulence. 

In this paper, we extend this work on turbulent sys- 
tems with ambipolar diffusion but without self-gravity. 
This paper has three main goals: First, we give a 
general discussion of the ambipolar diffusion Reynolds 
number, -RaDi that characterizes ambi polar diffusion 
(jMvers fc Khersonskvlll995l : IZweibelll2002f ) (§2). We then 
determine the numeric al values of Rab for the molecu- 
lar regions studied by iCrutcheJ (|1999D (§2) and show 
that they are consistent with the theoretically expected 
ones. Second, we use numerical simulations to determine 
the properties of the clumps that appear in a turbulent 
medium with ambipolar diffusion (§§3 & 4). In particu- 
lar, we show how the mass function and the mass-to-flux 
ratio of the clumps depend on Rab- Third, we ana- 
lyze the scaling properties of simulations with ambipolar 
diffusion and determine the range of physical parame- 
ters that characterize the simulations (§5). This discus- 
sion is continued in the Appendix, which gives a gen- 
eral discussion of how the results of simulations of tur- 
bulent boxes can be applied to physical systems, includ- 
ing those that are self-gravitating. Further results from 
these simulations, particularly those relevant to measur- 
ing the strength of the magnetic flcld and determining 
the effects of heating due to ambipolar diffusion, will be 
discussed in a future paper (Paper HI). 

2. THE AMBIPOLAR DIFFUSION REYNOLDS NUMBER 

The effects of ambipolar diffusion on a length scale i 
in a medium with a flow velocity v can be character- 
ized by the ambipolar diffusion (AD) Reynolds number, 
^ad(^)- This quantity appears to ha ve been first intro- 
duced bv lMvers fc Khersonskyl(|1995D : they referred to it 
as the magnetic Reynolds number, although that term is 
normally used to describe the effects of Ohmic resistiv- 
ity. The AD Reynolds number is motivated a s follows 
(IZweibel fc Brandenburd [19971: IZweibell [200l LMK): 
Ions in a partially ionized plasma are subject to two 
forces: the Lorentz force, ~ B'^^^^/A-kIb, where i?ims is 
the rms magnetic fleld strength and (.3 = |Srms/Vi?rms|; 
and the drag force, ^ab Pi PnV ad, where 7ad is the ion- 
neutral coupling coefficient, pi and p„ are the ion and 
neutral densities, respectively, and t;AD is the drift veloc- 
ity between the neutrals and the ions. When the ioniza- 
tion is low enough that the ion inertia can be neglected, 
these forces balance and the drift velocity is 



vad(^s) 



B^. 



^T^lABPiPn^B 



(1) 



We define the ambipolar-diffusion time over a length 
scale t as 



iAD(^) 



VAB{t) 



47r7ADP»p„^^ 
52 



(2) 



Similarly, we can introduce the ambipolar-diffusion 
length scale €ad, which is the length for which the am- 
bipolar drift velocity is the same as the flow velocity — 
i.e., in the frame of the ions, the length scale over which 
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the field varies in a steady flow: 



^AD 



Bl 



Dn collision tir 
the AD time scale and length scale are 



(3) 



In terms of the neutral-ion collision time, t^i — ^/^ADpi 



tAD - 



"^ pj"ni 



(4) 
(5) 



where ua = 5rms/(47rp)^/^ is the Alfven velocity and 
where we have assumed that the ion mass density is neg- 
ligible, so that pn ~ p. The efl^ect of ambipolar diffusion 
on a flow over a length scale (, with a characteristic ve- 
locity V is determined by the AD Reynolds number, 



i?AD(^) 






tf 






'i-^jADPiPn^V 



B^ 



(6) 



where tf = ijv is the flow time across a length i. Ob- 
serve that ambipolar diffusion increases in importance as 
-Rad(^) decreases; thus, it becomes more important at 
low densities, low ionizations, low velocities, small dis- 
tances and high held strengths. As [Myers & Khersonskv 
(|1995f l showed, the ratio of the size of the region, €, 
to the minimum wavelength of a propagating Alfven 
wave in which the inertia is provided b y the neutrals , 
Amin — T^vxtni (jKulsrud fc Pearcill969[ ). is directly pro- 
portional to RAr>{t}: 



A. 



ttMa ' 



(7) 



where A^a = v/va is the Alfven Mach number. 

We have defined the AD Reynolds number in terms of 
the mean densities and the rms field strength. In a su- 
personically turbulent medium, the densities are subject 
to large fluctuations, and if the Alfven Mach number is 
large also, the magnetic field has large fluctuations as 
well. If one deflnes Rad in terms of the local densities 
and field strength, then one can devise several different 
ways of averaging so as to obtain an effective value of 
i?AD for a turbulent medium; in particular, the length 
scale £ can be taken to be the size of the region or it can 
be determined self-consistently as the size of the aver- 
age eddy or density fluctuation. The resulting values for 
i?AD in a turbulent region of size £o range from slightly 
larger than i?AD(^o) to several times less (see Paper III 
for further discussion). One should thus bear in mind 
that i?AD(^o) is a characteristic value for the ratio of the 
ambipolar diffusion time to the flow time, and the actual 
value in a turbulent medium might differ from this by a 
factor of a few. 

Mouschovias (private communication) has emphasized 
that the AD Reynolds number is useful for turbulent 
media in which the velocity dispersion is determined by 
the turbulence (the case we are considering here), but 
not in systems in which the flow velocity is determined 
by the AD process itself. For example, in quasi-static, 
AD-regulated star formation, the AD length scale, -^ad 
is proportional to the radius of the self-gravitating cloud, 
and i?AD is of order unity. Similarly, Rab is not a useful 



parameter to characterize C-shocks (jDraind [19800 . since 
th e structure of such shocks adjusts itself so that Rad ^ 
1 (|Li. McKee. fc Klein|[200a . 

2.1. Numerical Evaluation oj Rad{£) 

Evaluation of the AD Reynolds number requires eval- 
uation of both the ion- neutral coupling coefficient, 7ad, 
and of the mean ionization mass fraction, Xi = Pil P- If 
this mass fraction is small {xi ^ l)i then pn c^ p and 



RABii) 



1 /tad/ 
2\ Cs 



X^M|3^ 



1ADP^\ XiM 



M 



(8) 

where c^ is the isothermal sound speed, M. = v/cs is 
the Mach number and /3 = Snpcl/Bf^^ is the plasma P 
parameter. We normalize our results to the case in which 
the ionization is dominated by HCO^. The ion- neutral 
coupling coefficient is then 



1.9 X 10-3 cm^ s-i 



7AD = 



m„ + mi 



^ 3.7 X 10" cm-^ s 



,3 ,.-1 „-l 



(9) 

(jDraine. Roberge. fc Dalgarnol [19831) . provided the rela- 
tive velocity of the ions and neutrals is less than about 19 
km s"^. Note that this value of the coupling coefficient 
differs from that adopted in LMKF due to our assump- 
tion that the ionization is dominated by IICO+. More 
generally, we shall write 



7AD = 3.7 X 10"7AD 



(10) 



where 7ad* is a number that allows for ions other than 
HCO+. 

Next, we consider the ionization. The processes that 
determine the ionization in molecular clouds are com- 
plex, and in general the ionization is time-dependent. 
We adopt a characteristic value of the ionization based 
on the assumption that the ionization is in a steady state 
and is dominated by HCO"*". In equilibrium, the mean 
ionization fraction by number is 



CcR 
aJT-H 



1/2 



(11) 



where njj = p/mh is the mean density of H nuclei, 
PYi = 2.34 X lO-^"* g is the mass per H nucleus, CcR ~ 
3 X 10^^^ s^^ is the cosmic ray ionization rate per 
H atom (see the discussion by iDalearnol l2006l) and a 
is th e relevant recombination rate (jMcKee fc Ostriked 
20071). Equat i on ([IT]) is consistent with the results of 
Padoan et al.l (j2004[ ) at late times and high densities 



for a = 2.5 x 10"^ cm'^ s"^, the value they adopted 
for the dissociative recombination rate of HCO+. If 
small PAHs dominate the ionization, then the disso- 
ciative recombination rat e is about 10 times smaller 
(jWakelam fc Herbstl l2008f) and the ionization s everal 
times larger. iTassis fc Mouschovias! ()2005l . I2007D . who 
included the effects of charged grains, adopted a disso- 
ciative recombination rate a = 1.0 x 10"® cm'^ s"^; their 
results for the ionization are approximately consistent 
with equation pT|) for densities nn ^ 10* cm~^. 

Inference of the ionization from observations gener- 
ally requires knowledge of the cosmic ray ionization 
rate and the density as inputs to the chemical models 
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used to interpret th e observations () Williams et al.llTQQSl : 
iPadoan et al.|[200^ . The ionization can be characterized 
by the parameter Ci defined by 



by 



= C,. 



CcR. 
7lH 



1/2 



(12) 



In the equilibrium model above, Ci = a ^'^, which 
is 630 cm'^/^ s^/^ for the fiducial case. iWilliams et alJ 
(|1998| ) found a median ionization Xe — 4.5 x 10~* (note 
that they normalized their results to H2 , whereas we are 
normalizing with respect to H). Their adopted ioniza- 
tion rate (CcR — 2.5 x 10^^^ s^^) and density (nn = 
5 X lO"* cm-3) correspond to fi c^ 2000 cm^/^ s^/^. 
More recently, IPadoan et al.l ()2004D have interpreted 
these data with time dependent models and infer lower 
values of the ionization and therefore Ci. They find 
Xe — 7.5 X 10~^ — 3.5 X 10~^ and attribute the higher 
values to the effect of FUV photoionization; for their 
assumed ionization rate (CcR = 6 x 10"^^ s~^) and den- 
sity (2 X 10^ cm~^), the implied value of the ionization 
parameter is Ci ~ 600 cm^/^ s^/^, fortuitously close to 
our fiducial va lue. The difference be tween the values of 
Cj_ inferred by IWilliams et al.l (|1998f ) and IPadoan et al.l 
((^04) is a reflection of the uncertainties that remain in 
determining the ionization in molecular clouds. 

To evaluate the AD Reynolds number, we require the 
ion mass fraction Xii which is related to the ion number 
fraction Xi = ni/hn by Xi = Ximi/fin — > 20.7xi, where 
the numerical evaluation is for HCO^. We then have 



Ci f CcR 



m 



n-H 



1/2 



2.25 X 10" 



^H.3/ 



where the numerical factor 



Q 



29 amu V 630 cm3/2 s^a / V 3 x lO-i^ s 



CcR 



(13) 



1/2 



(14) 

allows for deviations from the fiducial case and where 
fiH, 3 = nH/(10^ cm~3). Under the assumption that 
the ion mass is indee d about 29 amu (i.e., the mass of 
HCO+), the results of IWilliams et all (119 981 correspond 
to X* — 3, whereas those of (jPadoan et al. 2004 ) corre- 
spond to Xi — 1- The ionization can also be expressed 
in the form pi = Cp^/^ (Shu 1983), with C = XjP^^^ = 
1.09 X lO-i^x? g^^^ cm-3/2. Shu (1983) adopted a value 
C = 3 X 10^^^ g^^^, corres ponding to x* 7 ^ 3, in agree- 
ment with the estimate of IWilliams et all ()1998l ). Nu- 
merically, the AD Reynolds number is then 



RADii) = IG.OjAB* x:MI3 






T, 



1/2 



(15) 



where Ti = T/(10 K) and i^^ = (./{I pc). 



2.2. Regimes of Am,bipolar Diffusion 

We can distinguish several regimes in ambipolar diffu- 
sion in a turbulent medium. For Xi = Pi/ P "^ Ij we have 
p„ ~ p so that the neutral-ion collision time, t„i and the 
corresponding ion- neutral collision time, Un are related 



^in, — 



1 



Xi 



= XiU. 



(16) 



lADPn 1 AD Pi 

Similarly, the ion-Alfven Mach number, A^ai = 
v{4:npiy^'^ / Brras, is related to the Alfven Mach number. 
Ma, by Mai^ = Xi-^A^- It follows that 



Rad - Ma^ ( ^ 



MAi' ( j^ 



(17) 



We now identify five different regimes for ambipolar 
diffusion. For simplicity we ignore possible differences 
between the velocity dispersions of the neutrals and ions 
(to be discussed in Paper III), which could change the 
coefficient in front of Mai^ by up to a factor 2 in the 
expressions below. 

I. Ideal MHD {tf/tni -^ 00, corresponding to Rad ^■ 
00 for a given value of A^a): The ions and neutrals 
are perfectly coupled. 

II. Standard AD {tf > tni 3> tin, corresponding to 

-Rad > Ma^): The neutrals and ions are coupled 
together over a flow time so that the AD is weak. 
For M.A — C(l), linear Alfven waves can prop- 
agate, since the propagation co ndition for Alfven 
waves of wavelength I derived bv lKulsrud fc Pearcd 
fl969) is equivalent to i?AD(^) > ttTWa (eq. [7]). 
The wave damping is weak (i.e., Ttf < 1, where 
r = hk^v\tni is the damping rate for low-frequency 
waves) for i?AD(^) > 27r2. 

III. Strong AD {tni > tj > tin, corresponding to 

M.A^ > Rad > A^Ai^): The neutrals are no longer 
coupled to the ions in a flow time, but the ions re- 
main coupled to the neutrals. For A^a — C>(1), 
Alfven waves cannot propagate since Amin/^ > 
■k/Ma > 1 (eq. El. 

IV. Weakly coupled {tin > tf, corresponding to 
A^Ai^ > -Rad): The ions and neutrals are only 
weakly coupled and act almost independently. The 
damping rate for the high-frequency Alfven waves 
that propagate in the ions is F = l/(2fi„), so 
these waves are weakly damped in this regime: 

Ttf = -RAD(^)/(2A^Ai^) < 1 for waves of wave- 
length £. The Heavy-Ion Approximation is based 
on the assumption that the ion inertia is negligible 
and therefore does not apply to this regime (see 
below). 

V. Hydrodynamics {tf /tin ~> or Xi ~^ 0; corre- 
sponding to i?AD — >■ for a given value of A^a): 
The neutrals are not affected by the trace ions and 
act purely hydrodynamically. One can of course 
recover the hydrodynamic limit by letting i? — >■ 
so that AIa — >■ 00; in that case, i?AD is uncon- 
strained. It should be noted that the boundary 
between the hydrodynamic regime and the weakly 
coupled regime is a matter of choice; if one de- 
mands that the ions have at most a 1% effect on 
the neutrals, for example, then J?ad would have to 
be smaller than if one demands that the effects be 
limited to 10%. 
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It should be borne in mind that in all regimes except 
the last (where it is irrelevant), we have assumed that 
the ions are well-coupled to the magnetic field — i.e., the 
ion gyrofrequency is much larger than the ion-neutral 
collision frequency, flitin S> 1. Although we have defined 
the AD regimes for arbitrary values of the Alfven Mach 
number (provided B is large enough that flitm S> 1), 
this characterization of ambipolar diffusion is most useful 
when Ma ~ 0{1), as it generally is in molecular gas in 
the interstellar medium. 

2.3. Observed Values of Rab 

iCrutcheil ()1999D has summarized sensitive Zeeman 
measurements of magnetic field strengths together with 
other physical parameters, including the plasma /3 and 
the Mach numbers, for 27 molecular clouds. Of these, 
12 have only an upper limit on the line-of-sight magnetic 
field. Table 1 lis ts the values o f the parameters from 
Tables 1 and 2 in iCrutcheii ()1999:) that we use to com- 
pute the corresponding Rad using equation ([TS]) . We 
take the length scale io to be the cloud diameter. We 
use Crutcher's correction for projection effects on the 
magnetic field: Zeeman observations determine the line- 
of-sight component of the field, i?ios, and on average 
the value of B^ that enters the plasma-/3 parameter is 
3i?[Qg. We assume that the parameters describing the 
ion- neutral coupling and the ionization (7ad* and x*) 
are unity. From Table 1, we see that clouds with mea- 
sured field strengths have Rab ranging from a few to 
^ 70. Because the range of Rab is so large, we quote 
the logarithmic average, defined as 



{R 



AD, 



log 



-|^g(logflAD). 



(18) 



the logarithmic mean and dispersion of the AD Reynolds 
number in these clouds is (i?AD)iog = 17 ± 0.4 dex. 
Clouds that have only upper limits on the magnetic field 
have an average lower limit on the AD Reynolds num- 
ber of {Rad)io„ = 22; if we discard L889 as an outlier 
because of its unusually high Mach number {M — 7.3), 
the logarithmic mean is 18, which is comparable to that 
of the clouds with measured fields. We also include the 
Alfven Mach number, Ma, in Table 1. All the clouds 
have Rad > Ma^, implying that these clouds are in 
the standard AD regime (i j2.2p . The ratio of the size of 
the cloud to the minimum Alfven wavelength is in the 
range 2 — 15 for clouds with measured field strengths. 
We conclude that the effects of ambipolar diffusion must 
be considered in studies of molecular clouds, at least in 
those regions shielded from the interstellar radiation field 
so that X* = '^(l): in agr eement with studies extending 
back for many years fe.g.. lMouschoviaslll987| ). 

2.4. Predicted Rab and Implied Self- Gravity 

As we now show, it is possible to predict the AD 
Reynolds number for self-gravitating clouds that have 
an ionization of the form given in equation (jlSp. As a 
corollary, we show that the strength of self-gravity is not 
a free parameter in simulations of ambipolar diffusion in 
a turbulent medium. 

The importance of self-gravity in a cloud of radius Rq 
or in a simulation box of size io = 2Ro is determined by 



the virial parameter (jBertoldi fc McKeill992f) . 



GMo 2GMo ' 



where 



V3 



Mc, 



(19) 



(20) 



is the ID velocity dispersion in the cloud. The virial 
parameter is thus proportional to the ratio of kinetic to 
gravitational energy. We wish to treat both real clouds, 
which we approximate as effectively spherical, and tur- 
bulent boxes. Of course , real clouds are not spherical 
(jBertoldi fc McKe"ill992l give the generalization to ellip- 
tical clouds), but keeping track of these two cases pro- 
vides a gauge of the importance of geometric effects; fur- 
ther more, the spherical clo ud model has long been in use 
(e.g., iSolomon et al.l[l987| ). Let the area and volume of 
the cloud or box be 



*0i 



V = cvil 



(21) 
(22) 



where ca — (7r/4, 1) and cy — (tt/G, 1) for a spherical 
cloud and a box, respectively. The virial parameter then 
becomes 

Since pi = XiP, equation ([5]) for the AD Reynolds number 
can be rewritten as 



i?AD(4) 



(s;i-)""^»«-'"''^-'- '^^' 



which shows that the AD Reynolds number is determined 
by the ionization, the Alfven Mach number and the virial 
parameter. Insofar as the ionization is a function of the 
density, Rad{^o) will also depend on density. However, 
in the case of greatest interest, in which Xi ^ X*P~^, 
where x* is a number that is unity in the fiducial case 
(eq. [Ti)l . the AD Reynolds number is fixed at 



i?AD(4) = 19.77AD*x: 



6cyav 
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Ml 
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Molecular cloud cores and clumps with measured mag- 
netic fields are typically self-gravitating, with Ovir ~ 1, 
and have Ma ^ 1 (Crutcher 1999). As a result, for 
the fiducial values of 7ad and Xi, such regions have 
Rab{(-) '^ 20. This predicted value is in good agree- 
ment with the observed values discussed in i i2.3l which 
have a logarithmic mean of 17 and a dispersion of 0.4 
dex. GMCs as a whole have larger values of Rab{(-) 
since most of their mass is photoionized by UV radiation 
()McKeelll989[ ). so that they have a higher ionization than 
the cores and clumps within them (e.g., in an envelope 
of a GMC in which the ionization is dominated by 0+, 
the ionization parameter is x* ^ 10^). 

The importance of self-gravity in a magnetized medium 
can also be expressed in terms of the ratio of the mass to 
the magnetic critical mass, M$, which is the minimum 
mass that can undergo gravitational collapse. In terms 
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of the magnetic fiux, $ 
mass is 



BcA^O' the magnetic critical 



Af$ = c$ 



qTJ2 



(26) 



wher e c$ = l/27r for a cold sheet (jNakano fc Nakamural 
11978^ and w 0.12 for a cloud with a flux-to- mass dis- 
tribution corresponding t o a uniform field threading a 
uniform spherical clo ud (jMouschovias fc Spitzeri 119761 : 
iTomisaka et alJ[l988l ). For c$ = l/27r, the ratio of the 
mass to the magnetic critical mass is 



M*,o 



Mo 
Af$ 



1/2 



-> 
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Ma (27) 



(spherical cloudpS) 



which provides a simple relation between the Alfven 
Mach number, Ma and the two parameters describ- 
ing the importance of self gravity in a magnetized, tur- 
bulent cloud, ttvir and /i$.o- The ratio /i$.o is some- 
times written as the ratio of the observed mass-to- 
flux ratio to the critical o ne, (M/$)obs/(-Af/$)crit (e.g., 
iTroland fc CrutcheHl2008l ). Using equation dH]), we find 
that the AD Reynolds number is given in terms of /i$ o 
by 

i?AD(4) = 13.2 (^^ 7AD*xrXAA^*,o; (29) 

the factor in parentheses is unity for a spherical cloud. 
Gravitationally bound clouds that are both magnetized 
and turbulent have /i$^ o somewhat greater than unity 
since the gravity has to overcome both the turbulent mo- 
tions and the magnetic field (McKcc 1989). This expres- 
sion thus gives a similar result to that in equation (|25l) 
for Ma ~ 1 and fi^^ o — 1 — 2. 

These relations for i?AD(^o) can be inverted to give the 
values of the virial parameter and the ratio of the mass to 
the critical mass in terms of i?AD(^o) and Ma- In other 
words, a simulation of a turbulent box with ambipolar 
diffusion [which requires specification of i?AD(^o) and 
A^a] necessarily implies the strength self-gravity would 
have were it to be included: 
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7AD Xt 

Rad{Io[ 

i?AD(4 



TAD'Xr 



Ma' 

1 
M'a 



(30) 
(31) 



where we have set c^i = cy = 1, as is appropriate for a 
simulation box. For avir S> 1, the neglect of self-gravity 
is self-consistent. Parameter choices that lead to values 
of Ofvir *^ 1 and /^$,o ^ 1 are not self-consistent, since 
self-gravity w ould lead to turbulent motio ns that render 
"vir > 1 ('e.g.. lKlessen fc Hennebeilil2009f ). 

We emphasize that "implied self gravity" does not 
mean that simulations of ambipolar diffusion mimic the 
effects of self gravity. Rather, it means that the strength 
of the self-gravity, were it to be included, is not a free 
parameter provided the ionization parameter %* is spec- 
ifiedQ By contrast, a simulation of a turbulent box with 

^ Of course, the relation between a^ir and -Rad(^o) also depends 



ideal MHD is scale free; the density can be chosen so 
that self-gravity would be negligible if it were included. 
This freedom does not exist in simulations of ambipolar 
diffusion. 

3. SIMULATIONS 

In this paper, we extend the LMKF study of super- 
sonic turbulence with ambipolar diffusion, focusing on 
the physical properties of the clumps formed purely as 
the result of turbulent fragmentation with no gravity. 
LMKF performed a series of 256'^ simulations in a peri- 
odic box using the code ZEUS-MPAD to investigate tur- 
bulence statistics in non-ideal MHD without self-gravity. 
Like LMKF, we drove the turbulence with a fixed driv- 
ing pattern over the wavenumber range 1 < fc < 2 (where 
fc = fcphys^o/27r = io/^ is t he normalized w avenumber) 
using the recipe described in lMac Lowl (|1999H . The driv- 
ing maintained the 3D Mach number at A^ = 3, which is 
only mildly supersonic. The corresponding line-of-sight 
Mach number — i.e., the ID Mach number M/\/3 — is less 
than 2. The magnetic field was initially uniform, with 
a strength set by a plasma- /3 parameter of 0.1, corre- 
sponding to an Alfven Mach number Ma = 0.67; the 
turbulence is thus sub-Alfvenic. During the simulations, 
the volume-averaged magnetic field changed by less than 
10%, and as a result the volume-averaged value of /3 re- 
mained within 10% of its initial value. As shown in Table 
1, this value of /3 is close to the median of the 15 clouds 
with measured magnetic fields. 

The focus of our effort is to determine how the prop- 
erties of the clumps vary with i?AD(^o), so for now we 
discuss our results in dimensionless form; the physical 
conditions corresponding to these simulations will be dis- 
cussed in Sj5] below. We note, however, that for systems 
satisfying the linewidth-size relation, a Mach number of 
3 corresponds to a box size £o ~ 0.4 pc. Like LMKF, 
we considered values of Rau{£o) from 0.12, close to the 
hydrodynamic limit, to 1200, close to the ideal MHD 
limit. The run with i?AD(^o) = 12 has conditions similar 
to those in observed clouds; as we shall see below, if we 
assume that the simulated region satisfies the linewidth- 
size relation, its density would be nn ~ lO'* cm~'^. In 
i 34.2.H we show that the inertial range of the simulated 
turbulence extends over the range ^o/3— £o/20, so the AD 
length scale £ad = •^o/-Rad(^o) is in the inertial range for 
this run. For all the other runs, the AD length scale is 
outside the inertial range. The run with i?AD(-^o) = 1.2 
focuses on scales less than ^ad and has a lower value 
of the AD Reynolds number than any of the clouds ob- 
served byfCrutcher ( 1999) , most of which are gravitation- 
ally bound. If the simulation satisfied the linewidth-size 
relation, it would have a density of nn — 400 cm~'^, corre- 
sponding to an unbound cloud. On the other hand, the 
runs with i?AD(-^o) = 120, 1200 focus on scales greater 
than ^AD and have higher AD Reynolds numbers than 
any of the clouds with measured magnetic fields in that 
sample. The i?AD(-^o) = 1200 run could be applied to the 
outer parts of GMCs, where the ionization is dominated 
by C"*". The run with i?AD(^o) = 0.12 represents the 

on A^Ai but this is fixed in simulations of turbulent boxes with 
.^A ^ 1; by contrast, whereas observed clouds have definite values 
of X * ; it i s not necessary to specify this quantity in the simulation — 
see E31 
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transition to the hydrodynamic limit, and is primarily of 
theoretical rather than practical interest. 

All models were run for 3i/, where tj = ^q/Mcs is 
the flow time. In order to improve the statistics and the 
resolution, we re-ran the five models m3c2r-l [i?AD(^o) = 
0.12] to m3c2r3 [i?AD(^o) = 1200; see Table 2] in LMKF 
with the same initial conditions but using a 512^ grid. 
All the results reported in this paper are the result of 
simulations on such a grid. The total computing time for 
aU the models was - 600,000 CPU hours on the NCSA 
machine Abe using 512 processors. 

We made two principal approximations in our simula- 
tions. First, as discussed in the Introduction, we used 
the Heavy Ion Approximation (LMK), adopting an ion- 
ization xm = T^Xio and a corresponding ion-neutral cou- 
pling coefScient 7ad oc 7ad/''^, with TZ ~ 10^; here 
the tilde denotes quantities measured in code units (see 
ijA.3p . The key to the Heavy Ion Approximation is that 
even though each of these parameters differs from the 
actual value by a factor of 10^, the ion- neutral coupling 
is governed by the product of the parameters and has 
the correct physical value. According to the discussion 
in §2.2, the five AD models and the ideal MHD model 
span three regimes of AD as listed in Table 2, based on 
the initial i?AD(^o) in equilibrium. 

Our second principal approximation was in our treat- 
ment of the ionization. Simulations can be carried out 
with various assumptions about the ionization, includ- 
ing ion conservation, ionization equilibrium and time- 
dependent ionization. Following LMKF, we assumed 
that the number of ions is conserved, so that the value of 
Xi for the entire box is constant. The density is initially 
uniform, so that the initial ionization mass fraction, Xio i 
is the same everywhere; we took it to be 10"^. LMKF 
demonstrated that the results were the same as in the 
case of ionization equilibrium (basically because the time 
for a neutral to exchange momentum with an ion is small 
compared to the ionization time scale). More generally, 
the ionization is time-dependent. The ratio of the flow 
time, tf = IqI M.Csi to the characteristic ionization time. 



ii, 



i/CcR (see eq. [U), is large: 
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(32) 
-,(33) 



where a* = a/(2.5 x 10"^ cm^ s^^), Ccr = Ccr/(3 x 
10^^^ s""'^), and we have used equation (jTSl) . It follows 
that the molecular gas is typically very close to ionization 
equilibrium (although it is not necessarily close to chemi- 
cal equilibrium). In simulations, the relevant comparison 
is between the flow time across a cell, tf /Mg, where Mg 
is the number of grid cells in the length of the box, and 
the ionization time. For our runs, which typically have 
Ug = 512, we have i//(A/'giion,cq) — 6i?AD(^o) for fidu- 
cial values of the parameters. Ionization equilibrium is 
thus a good approximation for all the cases we consider 
except i?AD(4) = 0.12. 

To test our use of the approxi mation of ion conserva- 
tion (see also the Appendix in iLi et al.l 120081 ). we ran 
several 256^ models with time-dependent ionization for 
different values of Rad- We find that the properties of 



the clumps in these runs are within a few percent of those 
in the corresponding 256^ runs with ion conservation, 
with the exception of the ion density. In fact, the mean 
ion density in the entire box in the time-dependent case 
is less than that in the ion conservation case by up to a 
factor ^2. As a result, the value of the AD Reynolds 
number is reduced by a corresponding factor, as shown in 
Table 2. For large i?AD(^o), the gas is close to ionization 
equilibrium, so that pi ex p^/^. With this relation for the 
ion density, the mean ion density, and hence i?AD(^o)j 
are reduced by only a small amount compared to the 
case of ion conservation for the low Mach number we 
are considering if the density PDF is a lognormal with 
a wid th similar to that found by IPadoan fc Nordlrmdl 
(j2002[ ). For small i?AD(^o) the deviations from ioniza- 
tion equilibrium are larger, and correspondingly the dif- 
ference between the time-dependent and ion conservation 
results are larger as well. In this paper, however, we are 
exploring the effects of changing i?AD(^o) by orders of 
magnitude, so changes of < 2 do not affect our conclu- 
sions. 

4. PHYSICAL PROPERTIES OF CLUMPS 

The formation of high-density clumps is a natural 
outcome in simulations of highly supersonic turbulence, 
whether a magnetic field is included or not. Fur- 
thermore, high-res olution t urbulence simulations (e.g. 
ILi et a l. 2004; Pado an et al.| l2007') produce a mass spec- 
trum of clumps that qualitatively resembles the stel- 
lar initial mass function (IMF), with a peak at low 
mass and a power-law tail at high masses. Recent 
observations of rnolecular cores (e.g. [ Tachih ara et al.l 
[200llOnishi et al.ll2002tlAlves. Lombardi. fc L ada 200^ 
suggest a similarity between the stellar IMF and the 
core mass function. (We foh ow the terminology of 
iWilliams. Blitz, fc McKe3l2000l and use the term "core" 
to refer to the subset of clumps that are gravitation- 
ally boun d and will form a sta r or sm all multiple stellar 
system .) IPadoan fc Nordlrmdl ()2002D and IPadoan et al.l 
|2007|) have proposed a turbulent fragmentation theory 
for the IMF that relates the index of the velocity power 
spectrum to the slope of the higher-mass end of the 
clump mass spectrum. LMKF showed that ambipolar 
diffusion changes the velocity power index, and we con- 
firm that conclusion in Paper HI. If the turbulent frag- 
mentation theory is correct, we would expect a change 
in the slope of the higher-mass end of the clump mass 
distribution between the ideal MHD and the AD tur- 
bul ence simulations as welL (It should be noted that 
the iHennebelle fc ChabrieH[2008l theory leads to a much 
smaller predicted difference in the slope of the IMF in 
these two cases.) 

We use a CLUMP FIND algorithm, based on th e algo- 
rithm developed bv iWilUams. De Geus. fc Blitd (fl994) . 
to determine the clumps in our simulations. We de- 
fine "clumps" as connected regions with a density larger 
than the mean density of the turbulent box and will 
use the term "CIMF" for "clump mass function," reserv- 
ing "CMF" for "core mass function." This distinction 
is appropriate for our simulations since they do not in- 
clude self gravity. The density con t ours are separated by 
5p = 0.04p, which IPadoan etliLJ (|2007D found to work 
well in distinguishing distinct clumps. In order to infer 
the effects of AD on the CIMF, we require the clumps 
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to be resolved. As mentioned in LMK, ZEUS-MPAD 
needs at least 3 to 6 zones to accurately distinguish the 
effects of AD from those of numerical diffusion. There- 
fore, we require clumps to have at least 6 zones in the 
mean radius, unless otherwise specified; this requirement 
is validated in the resolution study of CIMF in §5.2.1. 
In implementing this resolution requirement, we define 
the effective radius as re = (3T4/47r)^/'^, where Vc is the 
volume of the clump is determined by summing the vol- 
umes of each cell in the clump that has a density above 
threshold; thus, for a porous clump, r^ is less than the 
projected radius of the clump (see §4.3.2|) . This approach 
to setting the resolution requirement eliminates small, 
very porous clumps, which have a lot of structure that 
is not well resolved. By varying 6p, we found that the 
number of clumps with mean radius larger than 6 cells 
does not change when the separation of the density con- 
tours is smaller than 4%, thereby justifying our choice 
of 6p. Before constructing the CIMF, we verify that the 
clumps defined in our simulations satisfy the heavy-ion 
approximation . 

4.1. The Heavy-Ion Approximation for Clumps 

The condition for the validity of the heavy-ion approxi- 
mation is Rad ^ A^Ai, where the ion Alfven Mach num- 
ber, A^Ai, is smaller than the total Alfven Mach number. 
Ma, by a factor (p^/p)^/^ < 1 (LMK). To calculate the 
AD Reynolds number of a clump, i?AD,c, we use the 3D 
density-weighted velocity dispersion of the neutral gas, 
•\/3ct„, inside a clump as the flow velocity and the mean 
diameter of the clump, dc — 2rc, as the length scale. The 
ion Alfven Mach number of the clump, A^Ai,c, is taken 
to be the rms value of Mai of all the cells in the clump. 
We can re-write the definition of Rab in equation (|6]) for 
clumps as 
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(34) 



In Figure [H we plot A^aj c versus Rad. c for models 
m3c2r-l, m3c2rl, and m3c2r3 at i = 3i/; the results for 
models m3c2r0 and m3c2r2 lie between the nearby mod- 
els. The data points all have Rad,c ^ AJai c even for 
model m3c2r-l, which has the smallest value of i?AD(^o)- 
We have verified that this is true at other times as well. 
LMKF found that the Heavy Ion Approximation was 
valid for a turbulent box provided RAD{Ki)/-M\i > 30, 
where iyi is the length scale for ion-velocity variations, 
which is generally significantly smaller than the size of 
the box. We do not know how iyi in the clumps com- 
pares with the clump diameters. If we assume that 
the two length scales are comparable, then the require- 
ment for the validity of the Heavy Ion Approximation is 
C'hia ^ 30. This is well satisfied for all the clumps ex- 
cept those in model m3c2r-l, which has i?AD(^) =0.12 
and is the most diffusive run. For this run, the box as a 
whole has Chia — 10, and the Heavy Ion Approximation 
is at best marginally satisfied. We have not observed any 
problems associated with this, however. 

Two interesting features of the results are worth not- 
ing. First, almost all the clumps have smaller values 
of Max^ than the box as a whole; this is expected be- 



cause of the linewidth-size relation. The few data points 
with slightly higher values of Mai^ are due to large sta- 
tistical fluctuation in the ion density in a few clumps. 
Second, we note that the distribution of the data points 
is roughly parallel to the power law Rad oc A^aj^ (the 
straight line). This is because the factor Chia depends 
on two quantities, the column density, PnDc, and the ve- 
locity dispersion ratio, Un/cff, each of which is almost 
independent of A^aj- 

4.2. Clump Mass Function (CIMF) 
4.2.1. Resolution: The Sonic Length and the Inertial Range 

In studying the properties of the clumps that arise 
in boxes with supersonic turbulence, two length scales 
are important: the sonic length, is, and the minimum 
scale for the inertial range, ^in,min, which corresponds 
to the wavenumber fcin^max = ^o/^in,min- The sonic 
length, which is defined by the condition that the rms 
turbulent velocity in a box of size is equal the sound 
speed, gives a characteristic scale for density fiu ctu- 
ations in a sup crsonic ally turbulen t medium (iPadoanI 
119951 : iVazquez^emadeni et al.l l2003f ) . The sonic length 
should be well resolved in numerical simulations since 
it is important to resolve these density fluctuations and 
the turbulent motions that produce them. The resolu- 
tion condition is Aa; ^ is, where Ax is the size of a 
grid cell; equivalently, in terms of the sonic wavenumber 
ks = in I is, we have i(j//S.x = Ng » ks- We assume 
that the turbulence in the box exhibits a linewidth-size 
relation of the fornfl 
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where id is the effective minimum driving scale; the corre- 
sponding wavenumber is kd = io/id- In our simulations, 
kd — 2, and we find that the average Mach number in 
boxes of size id is indeed very nearly equal to that for 
the entire box. At ~ 3. We also find (? ~ i for i?AD(^o) 
in the range 0.12-12; for i?AD(4) = 120, 1200, we find 



The sonic length in a simulation is then 
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Correspondingly, we have 
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for q = y and q 



respectively. This satisfies the 



resolution condition Afg ~ 512 ^ ks for _Rad(^o) < 12; 
for i?AD(^o) — 120, 1200, this resolution condition is 
only marginally satisfied. 

Before leaving the topic of the sonic length, we note 
that it can be inferred for actu al molecular clouds as well. 
For q= \ (the observed value- lHever fc Brun3l2004D , the 
sonic length is related to the linewidth-size parameter 

^ Note that lKrumholz fc McKei P005l ) defined the sonic length 
with respect to the ID turbulent velocity, (Jnt = Cs{£o/^s id)'i and 
adopted q = \', the two versions of the sonic length are related by 



3i/(2g)£^^ corresponding to l^ id = 3^a for q ■■ 
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CTpc (eq. [131) by 
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We define the inertial range of the turbulence as the 
range of wavenumbers over which the power spectrum 
is a power law in k. In our simulations, this extends 



> fc > 3, where k. 



G "'in, max 
3 



20 



over the range /cj, 

for our 512'^ simulations and ~ 10 for the 256"^ simula- 
tions reported in LMKF. For k > fcin,maxj numerical dis- 
sipation becomes increasingly important. Another way 
of expressing this condition is that with ZEUS, numer- 
ical dissipation becomes important at about 1/lOth the 



minimum wavenumber, k 



in, max 



0.1 X {Afg/2). It is 
desirable to have the sonic length in the inertial range 
ks < fcin. max, and this is satisfied for the 512^ simula- 
tions with -Rad(^o) < 12. Determining whether this con- 
dition is a general requirement for accurate simulations 
of supersonic turbulence is beyond the scope of this pa- 
per. We note that this condition becomes increasingly 
difficult to satisfy as the Mach number increases. 

Figure[2]shows the clump mass distribution for the case 
of i?AD(^o) = 1200 (close to ideal MHD), at resolutions 
of 512'^ and 256^. We can make an approximate relation 
between the clump masses and wavenumbers by associ- 
ating a wavenumber kc = ia/Dc, where Dc is the clump 
diameter, to each clump. The corresponding clump mass 
is approximately 
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where pc is the average clump density. For the high- 
resolution run, the mean density of the clumps within 
the inertial range (fee < 20) is pc — 2.6p (i.e., the mean 
density is 2.6 times the minimum clump density). The 
higher-mass part of the CIMF appears to be a power law 
(this is justified in ^4.2.2l below). Observe that the slope 
of the CIMF changes at log Mc ~ —4.3, corresponding 
to fee — 30 = 1.5fcin,max- In fact, the clumps with such 
a mass have Dc ^ 12 — 20 cells. This is similar to both 
the maximum wavenumber in the inertial range and to 
the sonic wavenumber, which are also shown in Figure [51 
to within a factor of 2. In order to determine whether 
either of these parameters is associated with the change 
in slope, we also plot the clump mass spectrum for the 
corresponding 256^ run, for which fcin,max = 10 (ver- 
tical dashed line) is reduced by a factor 2 whereas kg 
is unchanged. The results are clear: The break in the 
clump mass spectrum in the low-resolution run occurs at 
half the wavenumber as in the high-resolution one. Fur- 
thermore, there is no discernable effect associated with 
the sonic wavenumber, although it would be desirable to 
test this conjecture with both higher resolution simula- 
tions and for higher Mach numbers than A^ = 3, the 
value in the present simulations. It therefore appears 
that the dominant effect in determining the deviation of 
the ClMFfrom a power law is the numerical dissipation 
that sets in for wavenumbers k > kin max- The results 



of our simulations can therefore address only the higher- 
mass portion of the CIMF, with a minimum diameter of 
12 cells. Figure [3l shows the 3D spatial distribution of 
clumps, identified by CLUMPFIND with minimum di- 
ameter of 12 cells, from a snapshot of model m3c2rl. 

4.2.2. Implications for the Turbulent Fragmentation Model 
for the IMF 

As remarked above, the similarity between the core 
mass function and the stellar IMF suggests that the 
IMF may be defined during the formation of cores in- 
side mol ecular clouds. In the turbu lent fragmentation 
model of iPadoan fc Nordlundl ()2002[ ) and iPadoan et al.l 
(I2007D . the mass distribution of cores (i.e., gravitation- 
ally unstable clumps) has the form 
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where ax is the dispersion of the density PDF. The 
power- law index of the core mass function at high masses, 
r, is related to the index, n^,, of the velocity power spec- 
trum, P{k)dk ex k~"'''dk, by 



r 



(42) 



for the strong-field, ideal MHD case {B > Bd), and 



5 — 2n„ 



(43) 



for the weak-field, ideal MHD case (B < Bcr), wh ich in- 
cludes the hydrodynamic case (jPadoan et al.ll20Q7l ): here 
the critical magnetic field i?cr is defined by the condition 
that the postshock gas pressure be comparable to the 
postshock magnetic pressure. iHennebelle fc Chabried 
f|2008i ) have introduced an improved theory for the IMF, 
but we cannot com ment on the dif f erence s between their 
results and those of IPadoan et al.l (|2007f) since our sim- 
ulations do not include self-gravity. As noted above, the 
Hennebelle fc Chabrier theory predicts a much smaller 
difference in the slope of the IMF between the magnetic 
and non-magnetic cases than does the theory of Padoan 
et al. 

The simulations of lPadoan et al.l (|2007D do not include 
self-gravity. Based on the discussion at the beginning of 
fjsl below, we note that if one specifies the temperature 
and adopts a linewidth-size relation, then it is possible to 
fix one parameter arbitrarily. For the Mach number they 
adopted {M = 10), their box size of ^o = 6 pc is in good 
agreement with the linewidth size relation in equation 
(|53l) . However, they chose a density fin = 2 x 10^ cm~'^, 
which results in a virial parameter a^ir — 0.028, far lower 
than observed values. 

In contrast to the Padoan et al simulations, we 
have included an additional physical process — ambipolar 
diffusion — so that the strength of the self-gravity is deter- 
mined by the parameter governing that process, i?AD(^o)j 
as discussed in §2.41 For the value of A^a we have 
adopted (A^a = 0.67) and for the fiducial values of 
the ion-neutral coupling coefficient 7ad and the ioniza- 
tion parameter Xi^ 1^6 virial parameter of the box is 
ofvir = 41/i?AD(^o)^ (eq. [30l). This is unphysically low 
for i?AD(^o) ^6: in nature, large values of the AD 
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Reynolds number are accompanied either by large val- 
ue s of A^ A (whic h is unlikely according to the results 
of lCrutcheiill999l) or by larger ionizations than implied 
hy Xi = 1- The core Bonnor-Ebert mass — that is, the 
Bonnor-Ebert mass based on the turbulent pressure in 
the ambient medium (eq. IA50|) — in our simulations is 



M 

5.4 



-Rad(^o) 



0.84 



Mp 



"vir 1 



Ma 



(44) 



where the last expression is for fiducial values of the pa- 
rameters. Of the five AD models in Table 2, only model 
m3c2rl, with -Rad(^o) = 12 (comparable to the observed 
values), yields a physically plausible Bonnor-Ebert mass. 
We therefore do not attempt to put our clump mass func- 
tion in physical units here (physical units are discussed 
in fJS] below) . What we can study is the slope of the 
higher-mass portion of the CIMF, which is independent 
of the choice of units. 

Ambipolar diffusion could introduce two changes in the 
value of F: First, the relation between F and riy could 
be different, varying from the hydrodynamic relation to 
the MHD one as Rad{^o) goes from to cx). Second, 
as shown in LMKF, the value of n« also depends on 
-Rad(^o)- In the pure hydrodynamic case, Uy = 2 for su- 
personic turbulence, whereas in the MHD case the value 
of n^ is not precisely known and could depend on the 
plasma (3. 

In Figure 21 we show the ClMFs of models m3c2r-l, 
m3c2rl, and mSi. LMKF showed that the density cor- 
relation between data sets at different times approaches 
zero in a time slightly less than tf. Therefore, in order 
to build up the statistics, we use three data sets in each 
model run, at i ~ f/, 2tf and 3t/. Adding all the clumps 
together to form a single data set, we use reduced x^ fit- 
ting to determine the higher-mass slope of the core mass 
function, Egt. We face two problems in determining Ffit: 
First, we do not know the the range of masses to include 
in the "higher-mass" data, and second, we do not want 
our answer to depend on the size of the bins used in bin- 
ning the data. We begin by dividing the data into 20 
logarithmically spaced mass bins. To address the first 
problem, we carry out fits beginning with only the three 
highest-mass bins, and then steadily increase the num- 
ber of bins used in the fitting until the peak of CIMF is 
reached. Initially, the value of x^ drops as the number 
of bins increases, since more data are contributing to the 
determination of the slope. However, when the number 
of bins is large enough that the CIMF begins to deviate 
from a power law, the reduced x^ will increase. To ad- 
dress the second problem, we increase the number of bins 
from 20 to 40 in increments of 5 and adopt the value of 
Ffit with the smallest reduced x^ from all five sets of fit- 
ting. Usually, the slopes corresponding to the minimum 
reduced x^ from different total bin numbers are close to 
each other. The resulting slopes are listed in Table 3. 

In view of the noise fluctuations in the higher-mass 
range of the CIMF, we have performed a two-sample 
Kolmogorov-Smirnov (K-S) test to determine whether 
this part of the CIMF can be fit with a power law. The 
mass range extends from the highest mass bin to the 



breakpoint determined by the x^ fitting procedure de- 
scribed above. The null hypothesis is that the higher- 
mass end of the CIMF from the simulation has the same 
distribution as a power law. Our results show that the 
K-S test on all five AD models and the ideal MHD model 
fails to reject the null hypothesis at the 5% confidence 
level. The p- values of all the K-S tests with different bin- 
ning are between 0.49 and 0.97. We conclude that the 
higher-mass portion of the CIMF is statistically consis- 
tent with a power law. For the ideal mhd case (model 
m3i) and for -Rad(^o) = 12 (model m3c2rl), the power 
law extends over the entire inertial range. However, in 
the limit of low Rad{^o) (model m3c2r-l), the power law 
extends only over the upper half of the inertial range; 
higher resolution and/or more samples are needed to de- 
termine if the inertial range is consistent with a power 
law in this case. 

With these 512'^ models, the clump statistics are ade- 
quate to demonstrate that the higher-mass slopes depend 
on i?AD- If turbulent fragmentation is correct, this is no 
surprise because LMKF found that the spectral indexes 
of the velocity power spectra also depend on Rab ■ Here 
we draw on the results of Paper HI, which gives more ac- 
curate values of the spectral index for the velocity of the 
neutrals, n^„(fc), than LMKF (see Table 3). The trend 
of spectral index changing from an Iroshnikov-Kraichnan 
(Iroshnikov" 19631: lKraichnanlll965[ ) to Burgers spectrum 
([Burgers 1974) as one goes from large to small Rab is 
still clear, as reported in LMKF. 

In the limit of ideal MHD (model m3i), the higher- 
mass slope is Ffit — 1.21 ± 0.09 (see Figure 3), which 
agrees quite well with the prediction F = 1.18 from equa- 
tion (l42t with spectral i ndex riy = 1.45 ± 0.05. Note 
that iPadoan et al.l (|2007D get somewhat different results 
{uy = 1.9 and F = 1.4), but this is presumably due to the 
difference in fiow conditions: they have 7W = 10 and /3 = 
1, whereas we have Ai = 3 and = 0.1. In their hydro - 
dynamic simulations, iBallesteros-Paredes et al.l (|2OO60 
found that the shape of the CIMF depends on the Mach 
number of the turbulence, consistent with our result. If 
the shape of the CIMF is significantly affected by the fiow 
conditions, then the Padoan et al. (2007) model would 
imply that the IMF depends on the environment, since 
regions of star fo rmation do not all have similar physi- 
cal conditions. As IBallesteros-Paredes et al.l()2006[) point 
out, this could be problematic in view of observational 
support for an IMF that is approximately universal. 

For the model m3c2rl, which has Rad = 12, compara- 
ble to the observed value (§2.3), the higher-mass slope is 
Ffit = 1.43 ± 0.10, which is consistent with the Salpeter 
value. As noted above, in the limit of low -Rad (model 
m3c2r-l), we are unable to fit the data with a power 
law that extends over the entire inertial range; The slope 
for the high-mass portion of the range for which a fit is 
possible is Fat — 2.41 ± 0.14, which continues the trend 
that the slope increases as i?AD(^o) decreases. Since this 
slope applies to only part of the inertial range, however, 
we are unable to check the validity of equation (j331)j 
which relates the slope of the CIMF to the velocity power 
spectrur n in the weak field ca se. Consistent with the re- 
sults of IPadoan et al.l ()2007D . this slope is significantly 
greater than the Salpeter value of the higher-mass slope, 

r = i.35. 

Comparison of numerical simulations of turbu- 
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lence with either pure hydrodynamics or ideal MHD 
has shown tha t magnetic fields suppress fragmen- 
tation (e.g. iPassqt. Vazguez-Semadeni. fc Pouaiietl 



19951 iGammie et al j " 
Hennebelle fc Tevssieri 



[2003t iPadoan et ahl [2003 
120081 ) . We can see this ef- 
fect in our simulations by comparing the ClMFs of 
models m3c2r-l, in which the neutrals are almost purely 
hydrodynamic, and m3i, with ideal MHD. Figure |4] 
shows that the number of clumps in the low-i?AD models 
is greater than in the high-i?AD models, except at the 
higher-mass end: The total number of clumps with 
Dc > 12 cells from the three time snapshots in the 
quasi-hydrodynamic model m3c2r-l is 2093 (Table 3), 
whereas it is 1033 in model m3i. The total mass of 
clumps in model m3c2r-l is ~ O.IOGMq, whereas it is 
~ 0.092Afo in model m3i. On average, the mass per 
clump in the quasi-hydrodynamic model m3c2r-l is 
smaller than that in the ideal MHD model, which is also 
consistent with prior simulations. 

The turbulent fragmentation model for the IMF pre- 
dicts that the core mass function (CMF) (i.e., the mass 
function of gravitationally bound clumps) is the same as 
the clump mass function (CIMF) at high masses, and it 
is based on the assumption that the IMF is proportional 
to the c ore mass function (the la tt er is predicted in the 
work of iMatzner fc McKed [2000l ) . IPadoan et all pOOTi ) 
emphasize that the predicted CIMF for hydrodynamic 
turbulence is much steeper than for ideal MHD turbu- 
lence, and our work confirms this. Our work shows that 
there is a continuous variation in the higher-mass slope 
of the CIMF due to the effects of ambipolar diffusion, 
such that the fraction of stars born at high mass should 
increase with Rad- 

4.3. Mass-To-Flux Ratios 

Ambipolar diffusion plays an important role in the core 
collapse process when the clump mass is less than or com- 
parable to the magnetic critical mass (eq. E5)) . Observa- 
tionally, only a limited number of cores have measured 
mass-to-fiux ratios due to the difficulty in making pre- 
cise Zeeman measurements. Furthermore, observations 
give only the line-of-sight values for the magnetic field 
and column density, so the value of /x$ ex Af/$ for any 
particular core is nec essarily uncertain. From a study 
of 34 dark cloud cores. iTroland fc Crutche^ poM ) found 
an average value of/z$,c = 1.4 — 2.1 after allowance for 
projection effects; the smaller value is based on flattened 
clouds, whereas the larger one is for spherical ones. The 
median values are larger by about 20%. Observed cores 
are thus somewhat magnetically supercritical. 

4.3.1. Resolution Study 

In this section, we check the convergence of the mass- 
to-flux ratios of the clumps by comparing the results from 
the 256"^ and 512'^ runs for model m3c2r3. To carry out 
the resolution study, we consider only clumps that have 
a mass at least equal to the minimum mass of clumps 
with Tc > 12 cells in the 512'^ model; for the 256'^ run, 
this corresponds to re ^ 6 cells. For easy comparison 
of clump mass-to-flux ratios among models, we elimi- 
nate the dependence of the mass-to-fiux ratio on clump 
mass by plotting the ratio (Ai*,c/M*,o)/(-Mc/-Mo)^/'^ ver- 
sus Mc/Mq for both the 256'^ and 512'^ runs in Figure 



[SI Curve fitting shows that the slope of the 256^ data is 
0.06 ±0.03 and the slope of 512^ data is 0.07 ±0.02. The 
mean values of {ti<s>^c/ f'<!>,o) / {Mc/MqY^^ are 1.68 ± 0.03 
and 1.78±0.03 for the 256'^ and 512"^ models, respectively. 
We conclude that the mass-to-fiux ratios of clumps in the 
512"^ model are converged. 

4.3.2. Effect of R ad (io) on the Mass-to-Flux Ratio 

As discussed in §4.2.21 above, our choice of parame- 
ters allows us to study the effect of varying i?AD(^o) 
on the mass-to-flux ratios, but at the expense of con- 
sidering models that would be unphysical were grav- 
ity to be included: Equation (I3ip implies /i$,o — 
0.17i?AD(^o)/(7AD*X*) for our simulations, which is in 
the observed range only for the Rad (^o) = 12 case. What 
is of interest then is how the normalized values of the 
mass-to-flux ratio vary with Rad- For example, the ra- 
tio of the mass-to-flux ratio for an individual clump, ^it^^a 
to that for the entire box is 



M, 



tl"S>,c 



M*,o BcTtR] 



■c,± 



Mo 



Be) So So 



where Rc,jl is the radius of the clump normal to the fleld 
threading the clump, B^ and Sq is the mean surface 
density for the turbulent box. For the cases we consider, 
the mean field in the clump is close to the mean field 
of the whole box since the relatively small value of the 
Alfven Mach number, M.a — 0.67, leads to a relatively 
uniform field, as discussed in !j4.4l below. As a result, 
for most clumps the mass-to-flux ratios are just propor- 
tional to the surface densities. The ratio Sq/Sc is just 
the number of clumps along a flux tube. Furthermore, 
since the density of the clumps is typically a few times the 
threshold density (see below eq. HO]) and is thus approxi- 
mately constant, it follows that the mass-to- flux ratio in 
the clumps is proportional to the cube root of the clump 
mass: 

^0„c oc Sc oc pcRc (X Rccc M^^. (46) 

We now use our simulations to determine whether am- 
bipolar diffusion in a turbulent medium affects the mass- 
to-flux ratio in clumps, even in the absence of self-gravity. 
In order to ensure that the clumps we study are in the 
higher-mass, power-law regime of the clump mass distri- 
bution so that numerical effects are minimal, we choose 
a minimum clump mass that is above the threshold for 
the higher-mass regime in all cases. This minimum mass 
corresponds to a minimum clump radius of Tc = 6 cells. 

To determine how ambipolar diffusion affects 
the mass-to-flux ratio, we compute the value of 
(m*,c/m*,o)/(-^^cA^o)^^'^ for all clumps in each model 
and plot the results in Figure |6l The mean values of 
M*.c/m$.o for the three models are also tabulated in 
Table 3 and are shown as the horizontal lines in Figure 
[6l In all the models, the values of fj,^^c for the clumps are 
smaller than //^^o for the whole box due to fragmentation 
along flux tubes. This effect h as been observed in other 
MHD turbulence simulations jVazguez-Semadeni et al.l 
I2005at iTillevfc Pudnti [200I . The typical value of 
the normalized mass-to-flux ratio, /i*,c/M$.o ^ 0.1, 
is set by our resolution, since the number of clumps 

increases with decreasing size and iJ,g,^c scales as Mc ■ 
We observe from Figure [6] and Table 3 that (/x$,c)/a'*,o 
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shows a small systematic increase from the large Rab 
model to the small and moderate Rad models. 

This table also shows that the mean density of 
clumps, (pc), in the three models increases systemati- 
cally as Rad decreases. The dispersion in the values 
of /Li$,c//U$,o/(-Mc/Mo)^/'^ in Figure m shows a significant 
variation as .Rad decreases. The dispersions of mass-to- 
flux ratio (not mass-to-flux divided by M^'^) are given 
in Table 3. The dispersion for -Rad(^o) = 12 is almost 
twice that for i?AD(^o) — 1200. The larger dispersion of 
/z$_c and higher density of clumps at i?AD(^o) = 12 than 
at high i?AD(^o) suggest that material can more easily 
cross magnetic field lines as Rad decreases. A further 
decrease in -Rad(^o) to 0.12 results in a higher density, 
but increased fragmentation of the clumps reduces the 
dispersion somewhat. We conclude that, even in the ab- 
sence of self-gravity, ambipolar diffusion has an effect on 
the mass-to-flux ratios of clumps. 

4.4. Other Physical Properties of Clumps 

In this section, we summarize a number of other phys- 
ical properties of the clumps as functions of Rad.c 
in Figure [7] by comparing the two models m3c2r-l 
[Rad{^o) = 0.12, strong AD] and m3c2r3 [i?AD(4) = 
1200, strong ion-neutral coupling], which represent the 
two extremes of ion-neutral coupling among our simula- 
tions. Figure [7] gives side-by-side plots of the normal- 
ized clump radii, rc/io, the ion and neutral densities, 
Pi,c/po and Pn,c/po, magnetic energy density, Ub,c/Ub,o, 
clump mass, Mc/Mq, and ionization mass fraction, Xij 
for clumps for the two models. 

Figure[7^ shows that the normalized radii of the clumps 
in the strong-coupling model (m3c2r3) are, on average, 
larger than those for the strong AD model (m3c2r-l). 
This is a result of more fragmentation in the strong 
AD case. The largest radius in m3c2r3 is about dou- 
ble that in m3c2r-l. Since we require clumps to have 
a radius larger than 6 cells, there is a sharp trunca- 
tion in the size distribut i ons at j rV^n — 6/512 — 0.012. 
iVazquez-Semadeni et al.l ()2005b( ) found that clumps in 
the non-magnetic case were smaller than those in the 
ideal MHD case, consistent with our result. 

Note that the clumps in model m3c2r3 have a smaller 
range of -Rad, c because of the strong coupling between 
ions and neutrals. This is seen in all other properties as 
well. Figures Uh and [7h show the normalized mean ion 
and neutral densities of the clumps. The sharp bottom 
edge in Figures [7^ and [T]: is the result of the density 
threshold pc > po we chose in defining the clumps. The 
variations in ion density are much smaller in the strong 
coupling case than in the strong AD case. This is also 
reflected in the ionization mass fraction in Figure O- The 
ionization mass fraction of clumps in the strong coupling 
model is about constant, but that of clumps in the strong 
AD model varies by almost 3 orders of magnitude. In 
Figure [Tjl, the magnetic field is barely perturbed by the 
turbulence in model m3c2r-l because of weak coupling; 
the magnetic field energy density in the clumps, Ub,c, is 
very nearly the same as that for the whole box. Although 
the magnetic field is perturbed more in model m3c2r3, 
most clumps have Ub,c within 50% of that in the box. 
Figure [Tjs shows that the largest clumps in m3c2r3 are 
more massive than the largest ones in m3c2r-l. This 



larger mass is due to a larger size, since the densities in 
the two models are about the same, and can be under- 
stood as the result of magnetic suppression of fragmenta- 
tion, as discussed in §5.2.2. The clump properties shown 
in Figure [7] include clumps down to re = 6 cells. From 
this figure, we see that the global physical properties of 
clumps scale smoothly from re = 6 cells to the largest 
clump. 

5. PHYSICAL UNITS FOR SIMULATIONS OF TURBULENT 
BOXES WITH AMBIPOLAR DIFFUSION 

The results of our simulations have been reported in 
dimensionless form. How can they be converted to 
physical values? A simulation of an isothermal, mag- 
netized, turbulent box is characterized by three dimen- 
sional parameters — the size of the box, £q, the mean den- 
sity in the box, riu, and the isothermal sound speed, 
Cg — and two dimensionless ones — the 3D sonic Mach 
number, M — 3^/^(T nt/cs and the plasma- /J parame- 
ter, 13 = Sirpcl/B'i (lOstriker. Gammie. fc Stond 119991: 
iPadoan fc Nordlundlll999[ ). Here nw is the mean densitv 
of hydrogen nuclei, ant is the ID nonthermal velocity 

1 In 

dispersion and Brms = {B"^) is the rms magnetic field. 
In the absence of other physical processes, all these pa- 
rameters can be selected arbitrarily, although the value 
of Cs (x T^/^ is tightly constrained for molecular clouds, 
which generally have temperatures in the range 10—20 K. 

Inclusion of a new physical process, such as ambipolar 
diffusion, introduces a new dimensional constant, in this 
case the ion- neutral coupling parameter, 7ad- Corre- 
spondingly, a new dimensionless parameter (in this case, 
Rab) can be formed and the number of independent di- 
mensional parameters is reduced by one. For a given 
sound speed, there is thus one independent dimensional 
parameter, such as the density, in simulations of ambipo- 
lar diffusion; such simulations are therefore scale free. 
Treatments of ambipolar diffusion require specification of 
the ionization, which in principle can introduce another 
dimensionless parameter that in tu rn w ould determine 
the scale. However, as discussed in ^A.3\ the Heavy Ion 
Approximation eliminates this constraint. 

Adoption of a linewidth-size relation, as is observed in 
molecular clouds (Larson 1981), also reduces the number 
of independent dimensional parameters by one. Hence, 
if an isothermal system satisfies a linewidth-size relation 
and is subject to ambipolar diffusion, then its velocity 
scale is set by the isothermal sound speed, Cg oc T^/^, 
and its size and mean density are determined by dimen- 
sionless parameters. In this case, a given simulation ap- 
plies to only one set of parameters describing the box. 
This is discussed further in the Appendix, which gives 
explicit expressions for properties of turbulent boxes in 
the general case, when they satisfy a linewidth-size re- 
lation, and for self-gravitating boxes. Here we present 
the scaling for our simulations of turbulent boxes with 
ambipolar diffusion. 

5.1. General Scaling Relations 

To determine how simulations of a turbulent box with 
ambipolar diffusion can be scaled to physical systems, we 
use equation ([T5|) to solve for the size of the simulation 
box, Iq. We find that it is determined by the remaining 
two dimensional parameters (fin and T) along with five 
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dimensionless parameters [i?AD(^o)i 



£o = 0.031 



7AD*X» 



M f Ti 



7AD , Xi 1 
1/2 



Ma^ V"-h. 



13 and M]: 
(47) 



pc, 



when 

Nn < iVLWS = 1 
or, equivalently, when 



3 X lO^V* 



The flow time across the box, the mass in the box, and 
the column density are then 



"H < n-LWS = 9.6 X 10 



'pc 



tf 



Mcs 



Mo = p£\ 



1.62 X 10^ 



1.06 X 10' 



Rab{ 



'Q 



Tad'X, 
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^H = nH4 = 9.6 X 10 
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For larger values of the column density and density, 
the linewidth-size relation must take the effects of self- 
gravity into account. The resulting virialized linewidth- 
\^^)size relation has a ex (Si?)^/^ and is equivale nt to settin; 



M/ 



, fi 

the v irial parameter equal to unity, avir = 1 (IHever et al.l 

^,^ 1 20091 : see §A.2.1| . The linewidth is greater than that in 

('T-H,3-'i) cn^SaJ^j^g turbulence-dominated case due to the effects of self 



Note that these scalings are preserved by the Heavy-Ion 
Approximation, in which the ion mass fraction (ex x*) 
is increased and the ion-neutral coupling coefficient (ex 
7AD*) is decreased by the same factor. The strength of 
the magnetic field does not depend on i?AD(^o): 



B = {Awpc 



.2U/2 



M 



3.2{nn,3T, 



,1/2 



M 



HG. (51) 



Ma Ma 

As discussed in i j2.4[ simulations of gas in which the 

— 1/2 

ionization scales as rijj have an implicit value of the 
virial parameter, avir, given by equation pop . Actual 
physical systems have avir ^ 1, since violations of this 
inequality lead to gravitational motions that raise avir 
up to order unity. Hence, this sets a lower limit on the 
product of the ionization parameter and the AD coupling 
parameter for a given value of i?AD(^o)j 



7AD X, 



> 



i?AD(4) 



(52) 
14.2Xa^ 

The lower limit on the ionization corresponds to the case 
of gravitationally bound clouds and clumps discussed in 
(for spherical clouds, the coefficient 14.2 is replaced 

19.7). 



by 14.2/ct/' 



5.2. Scaling with the Linewidth- Size Relation 

Most molecular gas in the Galaxy is observed to obey 
a linewidth-size relation 



fnt 






(53) 



where Rpc is the radius of the region measured in pc 
and typically Cpc — 0.72 km s~^ (jMcKee fc Ostrikeil 
l2007f) . The linewidth-size relation is quite general: it 
applies to within a factor ~ 3 to molecular gas ranging 
from small clump s much less than 1 pc in size to GMCs 
(jFalgarone et al.ll2009D . Taking R = £o/2 and noting 
that M is the 3D Mach number, we find 
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ol/2 ^ 
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where 



'pc 



_ '^pc 

" 0.72 km s-^ 
([20I0I) have 



iFalgarone fc McKei 
turbulence- dominated linewidth-size relation applies only 



(54) 
(55) 
shown that this 



gravity. 

When the turbulence-dominated linewidth-size rela- 
tion applies, so that TVh and nn satisfy the inequalities in 
equations (|56|) and ()57p . then the size of the simulation 
box is determined by equation (j54p : 



£0 = 0.0454 



M^Ti 



pc. 



(58) 



With the aid of equation (|47)) . one can then express the 
density in terms of the linewidth-size parameter, a*, 



nH = 9.6x 10* 
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RAuilo) 



14.27Wa'7ad* 
2 



14.2Xa'7ad*X* 



(60) 



where the factor in brackets is < 1 since the correspond- 
ing virial parameter must be > 1 (eq. I52p . Similarly, one 
can show that 



N» 



N, 



LWS 



RAuiio) 



U.2Ma'iab*X* 



(61) 



with the aid of equation (j56p . The mass corresponding 
to iVLws a-iid ^Lws — i-e., the maximum mass at which 
the turbulence-dominated linewidth s ize re lation holds — 
is Mlvvs I which is given in equation (JA41P . The mass in 
the simulation box is given in terms of Mlws by 



Mo = Mlws 



RAuiio) 



U.2Ma'iad*X* 



(62) 



On the other hand, when the system being simulated 
is self-gravitating, then avir ~ 1 and the inequality in 
equation ([5^ is replaced by an equality. Equations 



(jiTj) . (|49l) . and (EOl) show that in this case, £, 
and Mo oc n 



'1/2 



oc njj 



1/2 



'H 



are smaller than in the turbulence- 



-1/2 . 



dominated case, whereas A'h oc n^ is larger. The gen- 
eral case is discussed in the Appendix, !)A.2.3I 

5.3. Physical Parameters for Simulations 

We are now in a position to discuss the physical param- 
eters corresponding to our simulations. For simplicity, 
we shall assume that the temperature is T = 10 K and 
that the linewidth-size parameter has its fiducial value. 
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CTpj, — 1, corresponding to cTpc = 0.72 km s^^. The 
maximum column density for the turbulence-dominated 
linewidth-size relation is -/Vlws = 1-3 x 10^^ cm~^, and 
since the Mach number is A^ = 3, the corresponding 
maximum density is riLws = 1-1 x 10''' cm"''. We shall 
focus on the four cases i?AD(^o) = 1-2, 12, 120, 1200, 
since the i?AD(^o) = 0.12 case was done to model the 
transitio n to the h ydrodynamic limit. Recall that the 
clouds in lCrutcheii (1999)'s observations have a logarith- 
mic mean value (i?AD)io(>; = 17, comparable to the value 
in the i?AD(^o) = 12 simulation. The RAui^o) = 1-2 
simulation has a somewhat smaller value of -Rad, and 
the i?AD(^o) — 120 simulation a somewhat larger value, 
than any of the clouds in that sample; however, it must 
be borne in mind that this sample by no means covers all 
the types of molecular gas in the Galaxy. In particular, 
the i?AD(^o) — 1200 run is relevant to the outer parts of 
GMCs, where the ionization is dominated by C"*". 

The virial parameter associated with a given value of 
^ad(^o) in our simulations is 



RabJJq) 
6.47AD*xr 



(63) 



from equation (|30p . Since our simulations have A^a = 
0.67, the constraint on the ionization set by the require- 
ment ttvir ^ 1 implies (eq. I52p 

(64) 



RAuiio) < 6.47AD*X*- 



First consider the case in which i?AD(^o) = 1-2. 
We assume that the ionization and coupling param- 
eters have their fiducial values (7ad* — X* = 1); 
the ionization constraint is then well satisfied. The 
virial parameter is avir = 28 from equation ()63|) . so 
the self-gravity is negligible in the system being simu- 
lated. Equations (gT]), dH]), ^ and dH]) imply that 

- —1/2 

the size of the system is ^o — 0.25njj 3 pc, the mass 

-—1/2 

is M = 0.54rijj 3 M©, the column density is A^h — 

7.7 X lO^^ri-jj 3 cm~^, and the magnetic field is i? = 

- 1/2 
14njj' 3 /iG. Much of the unbound molecular gas in the 

Galaxy s atisfies the turbulence -dominated linewidth-size 
relation (jFalgarone et al.|[2009() . If the simulated system 
satisfies this relation, then the density is fin = 370 cm"'' 
from equation (I60p , and correspondingly the size of the 
simulation box is £0 = 0.4 pc, the mass is Mq = 0.9 Mq, 
the column is A'h = 4.7 x 10^" cm~^, and the magnetic 
field is S = 8.7 /iG. 

Next, consider the simulations with i?AD(^o) = 
12, 120. For these runs, the implied virial parameter 
is less than unity for the fiducial values of the ioniza- 
tion and coupling parameters. If 7AD*Xi' is as close as 
possible to its fiducial value, then the inequality in equa- 
tion ([M|) becomes an equality, and equations (|T7)) . (gH]), 

and imply 4 = l-3fiH,^3^ PC, Mq = ^2n^]{^ Mq, 
and A^H = 4.1 x lO^^rijj 3 cm"^. These conditions cor- 
respond to avir = 1, so the systems are on the virialized 
linewidth-size relation. As remarked above, the virialized 
linewidth-size relation applies when the density and col- 
umn density are large, tj-h ^ '^lws = 1.1 x 10'^ cm"'' and 
^H > ^Lws = 1.3 X 10^^ cm~^. Correspondingly, the 



size of the system is €0 ^ 0.4 pc the mass is Mq < 
-^LWS = 25M0, and the magnetic field is i? > 50 /xG. 
In sum, the systems we have simulated are relatively 

1 /9 

small, with £q < l/n^ 3 pc (eqs. ST] and [SU and 

Mq < lOO/nJ/3 Mq (eq. mi). If the system being simu- 
lated lies on the linewidth-size relation, then its mass is 
M < Mums = 25 Mq. As shown in gA.2.3[ this inequal- 
ity also holds if the system has a linewidth greater than 
that given by the linewidth-size relation. Similarly, equa- 
tion (|47)) shows that the size of the system decreases with 



cr* ; as a result, if the system being simulated lies on or 
above the linewidth-size relation, then its size is no larger 
than the size corresponding to the turbulence-dominated 
linewidth-size relation, £0 < 0.4 pc. Reference to equa- 
tions (|Tf)) and (|49l) shows that simulations of larger re- 
gions would require higher Mach numbers, given that the 
Alfven Mach number is observed to be of order unity. 

6. CONCLUSIONS 

Ambipolar diffusion is a key process in molecular 
clouds since it redistributes magnetic fiux and damps 
waves. The importance of ambipolar diffusion in a tur- 
bulent medium on a length scale £ and velocity disperion 

V is governed by the AD Reynolds number Rab — £/£a'd, 
where £ad — v\tni/v is the length scale over which the 
magnetic field must vary in orde r to have a drift velocity 

V between the neutra ls and ions ( Zweibel fc Brandenburg 
119971: I Zweibel [2001 . [Note that Rad is useful in de- 
scribing ambipolar diffusion whenever the velocity field 
includes a significant turbulent component; it is not use- 
ful for non-turbulen t, AD-driven gravitational collapse 
(jMouschoviasI I1987D or C-shocks (jDraind [T980l ) . where 
Rad is of order unity.] We have carried out two-fluid 
simulations of isothermal, turbulent boxes using the code 
ZEUS MPAD (described in LMK) at a resolution of 512^ 
for AD Reynolds numbers ranging from -Rad = 0.12 to 
i?AD = 1200, plus a simulation with ideal MHD. The 
resolution we have used is suflicient to resolve the sonic 
length within the inertial range, permitting accurate sim- 
ulations for our A^ = 3 calculations. The mean Mach 
numbers were fixed at A^ = 3 and A^a = 0.67, corre- 
sponding to a plasma- /3 parameter /3 — 0.1. The pur- 
pose of our simulations was to determine how the prop- 
erties of the clumps formed in molecular clouds depend 
on i?AD(^o)- One of the simulations (with i?AD(^o) — 12) 
was in the middle of the observed range of the observed 
values of the AD Reynolds number; two of the simula- 
tions (those with -Rad(^o) — 1-2 and 120) were somewhat 
below and above the observed values of i?AD(^o); and 
the remaining two simulations, with i?AD(^o) = 0.12 and 
1200, were designed to show the transition to hydrody- 
namics and ideal MHD, respectively. In order to carry 
out these simulations, we used the Heavy Ion Approxi- 
mation with an ionized mass fraction of 10"^ (LMK) to 
represent physical systems with actual ionized mass frac- 
tions ~ 10"^. We validated our simulations with conver- 
gence studies at lower resolution. The power spectra in 
our simulations show that the inertial range of our sim- 
ulations extends down to a length scale £q/2Q, which is 
comparable to the sonic length; it is important to resolve 
the sonic length in simulations of turbulent boxes since 
the density has significant fluctuations on larger scales. 
Our principal conclusions are: 
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1 . Values of the AD Reynolds number Rau in a sam- 
ple of 15 n iolecular clump s with measured mag- 
netic fields ()Crutcheiill999l ) range from 3 to 73; the 
logarithmic mean value is 17. Omitting one out- 
lier, the clumps with upper limits on the magnetic 
field have an average lower limit of (i?AD)iog > 18- 
The predicted value of the AD Reynolds number for 
self-gravitating molecular clouds and clumps with 
the fiducial ionization is Rad — 20, in excellent 
agreement with observation. 

2. Several regimes of ambipolar diffusion can be iden- 
tified, depending on the ratio of the flow time, tf, to 
the ion- neutral collision time, tin, and the neutral- 
ion collision time, tni'- (I) ideal MHD {tf/tni —>■ oo, 
corresponding to Rau — > oo for a given value 
of A^a); (II) standard ambipolar diffusion, with 
tf > tni, SO that the neutrals and ions are cou- 
pled together over a flow time; (III) strong AD 
(tni > tf > tin), so that the neutrals are not 
coupled to the ions over a flow time, but the ions 
are coupled to the neutrals; (IV) weakly coupled 
{tin > tf), SO that the ions and neutrals behave al- 
most independently over a flow time; and (V) hy- 
drodynamics {tf /tin — >■ or Xi — >• 0, corresponding 
to Rati — >■ 0). The molecular clumps in Crutcher's 
sample are all in the second regime, standard AD. 

3. Implied self-gravity: Since the ionization scales 
approximately as the square root of the density, 
the ambipolar diffusion ti me is proportional to the 
gravitational free-fall time (j Mouschoviasll 1 9 87l) . As 
a result, any simulation of ambipolar diffusion has 
a gravitational virial parameter ctvir that is de- 
termined by Rau, -^a and the parameters de- 
scribing the ion-neutral coupling and the ionization 

[ofvir oc (7AD*X*/i?AD)^A^A''— eq. [SOI- It is not 
possible to carry out a simulation in which the ef- 
fects of self-gravity and ambipolar diffusion arc var- 
ied independently unless the ionization is treated as 
a free parameter. 

4. Clump mass spectrum. Us ing Clumpfind 
()Williams. De Geus. fc Bhtg IT991 . we found all 
the clumps with densities exceeding the mean den- 
sity in the box. We find that the slope of the higher- 
mass portion of the resulting clump mass spec- 
trum increases as Rau decreases, which is quali- 
tatively consistent with Padoan et al's (2007) find- 
ing that the mass spectrum in hydrodynamic tur- 
bulence is significantly steeper than in ideal MHD 
turbulence. The value of the slope that we find for 
Rau = 12, the case closest to the value observed 
in molecular clouds, is Ffit — 1.43 ± 0.10, which 
is consistent with the Salpeter value, T = 1.35. 
The almost-ideal MHD case {Rau = 1200) has 
a slope Fgt = 1.22 ± 0.11, which is marginally 
consistent with the Salpeter value. We further 
confirm Padoan et al's (2007) relation between 
the index of the power spectrum and the slope 
of the clump mass spectrum in the limiting cases 
of ideal MHD and near hydrodynamics. How- 
ever, the value we find for the spectral index in 
our ideal MHD simulation differs from theirs, pre- 
sumably because our simulation has lower values 



of 13 and M . This s uggest s that the IMF in the 
iPadoan fc Nordlundl pOOl : iPadoan et all pOOTi ) 
turbulent fragmentation model depends on the en- 
vironment, which could conflict with evidence for 
an IMF that is approx imately universal (see also 
IBallesteros-Paredes et al. 2006) . 

5. Ambipolar diffusion affects the mass-to- flux ratio 
of clumps, even in the absence of self- gravity: The 
average mass-to-flux ratio /i$.c/A'*,o at low Rau is 
slightly larger than at high Rau , and the dispersion 
in the values of /i$.c/M*,o for individual clumps at 
moderate -Rad is almost twice that at high i?AD- 

6. Scaling relations for simulations of isothermal tur- 
bulent boxes. A simulation of an isothermal, mag- 
netized, turbulent box is characterized by three 
dimensional parameters: the size of the box, £oj 
the mean density in the box, p oc tt-h, and the 
sound speed, c^ (e.g., lOstriker. Gammie. fc Stong 
[1999). A single simulation with ideal MHD ap- 
plies to an infinite range of values of each of these 
dimensional parameters, provided that the dimen- 
sionless parameters describing the simulation (in 
this case, t he Mach numbers .M and A ^a) are the 
same (e.g.. IPadoan fc Nordlundl I1999D . Except in 
regions of high-mass star formation, molecular gas 
generally has a temperature T ~ 10 — 20 K, so that 
Cs is nearly constant; as a result, there are only 
two dimensional parameters that have a significant 
variation, (.q and p. Each physical process that 
is introduced into the simulation, such as ambipo- 
lar diffusion, introduces a dimensionless parameter, 
such as Rau, which must be fixed for the simula- 
tion, thereby reducing the number of scaling pa- 
rameters by one. For simulations with ambipolar 
diffusion, the physical parameters describing the 
system being simulated are characterized by a sin- 
gle dimensional parameter (for constant Cs), which 
we took to be the mean density (^STT]). Such sim- 
ulations are intrinsically scale free. Even if one 
includes self-gravity, then, so long as the ionization 

scales as Xi oc n^ , the simulation remains scale 
free. However, if one further requires that the sim- 
ulation satisfy an observed linewidth-size relation, 
then the mean density is determined and there are 
no independent scaling parameters. 

7. Physical parameters associated with the simula- 
tions. Two of the simulations we carried out were 
for the purpose of studying the transition to ideal 
MHD [Rau{£o) = 1200] and to hydrodynamics 
[i?AD(^o) = 0.12], so we focus on the remaining 
three, with -Rad (4) = 1-2, 12, 120. As discussed 
in item (3) above, a simulation of ambipolar diffu- 
sion has an associated value of the virial parame- 
ter. The virial parameter cannot be significantly 
less than unity since the self-gravity would induce 
collapse that would lead to Ofvir ^ 1. This con- 
straint gives a lower bound on the ionization such 
that -Rad(^o) ^ 6.47AD*X*- As a result, two of 
the simulations [RAu{io) = 12, 120] could be re- 
alized in nature only if the ionization and/or the 
ion-neutral coupling constant were larger than the 
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fiducial values, which correspond to 7AD*Xi = 1- 
The -Rad(^o) = 1-2 simulation corresponds to a 
system in which self-gravity is unimportant (for 
7AD*Xi* = l)i whereas the i?AD(4) = 12, 120 sim- 
ulations correspond to systems that mostly likely 
are gravitationally bound (for 7ad*X* ^^ close to 
unity as possible). Because of this constraint and 
because of the small value of the Mach number we 
adopted {A4 =3), our simulations apply to small 

1/2 

regions in molecular clouds, with £o < l/n-^ 3 pc 

1 /9 

and A/o < lOO/fijj' 3 A/q. If the system being sim- 
ulated has a velocity dispersion on or above the 
linewidth-size relation observed in the Galaxy, then 
the size of the region is £0 < 0.4 pc and the mass 
is < 25 Mq. 



A general discussion of scaling relations for self- 
gravitating systems is given in the Appendix. 
In applying the linewidth-size relation, we follow 



iFalgarone fc McKed ()2010f ) in distinguishing the 
turbulence-dominated relation from the virialized 



one. 
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APPENDIX 
SCALING LAWS FOR ISOTHERMAL TURBULENT BOXES 

In this Appendix, we give a general discussion of scaling laws for simulations of isothermal, turbulent gases in a box. 
Although we do not include the effects of self-gravity in the text, we do include it here, so as to make the discussion 



more generally useful. We focus on molecular gas, since such gas is generally approximately isothermal. Particular 
scaling relations that have been derived previously are not ed 
iTiUev fc Pudnt3[200l I Vazquez-Semadeni et al.ll2005aL I2OO8). 



Imml 



Ostriker. Gammie. fc Sto^ 119991: 



In the simplest case in which there is no gravity and the M HD is ideal, a simulation o f an isothermal, magnetized, 
turbulent box is characterized by two dimensionless parameters (jPadoan fc Nordlundlll999l ). the 3D sonic Mach number, 
A4 = 3^/^o-nt/cs, and the plasma-/3 parameter, /3 = S-Kpc^/Bf^^. Here a^t is the ID nonthermal velocity dispersion, Cg 

1 /2 

the isothermal sound speed, p the mean mass density and Sims = (-B^) the rms magnetic field. Equivalently, the 
two parameters can be chosen to be the sonic Mach number and the Alfven Mach number, A^a = 3^/^crnt/wA, since 
the plasma- /3 parameter is related to A4 and A^a by 



13 = 2 



Ma 
M 



(Al) 



In general, the Mach numbers, A^ and A^a, and the plasma- /3 parameter are functions of time. 

The turbulent box is also characterized by three dimensional parameters: the size of the box, (,q, the mean 
density in the box, p = A /p/^n, where Mq is the mass in the box, and the isothermal sound speed, Cs (e.g., 
lOstriker. Gammie. fc Stoiiail999) . In the absence of other physical processes, these parameters can be selected arbi- 
trarily. In other words, a given simulation corresponds to definite values of M and /?, but it can be scaled to arbitrary 
values of p, £o: and Cg. However, the introduction of a new physical process, such as self-gravity or ambipolar diffu- 
sion, introduces a new dimensional constant and a corresponding new dimensionless parameter, so that the number of 
independent dimensional parameters is reduced by one. The same reduction occurs if a relation between dimensional 
parameters is assumed, such as a relation between the size of the box and the mean velocity dispersion (a linewidth- 
size relation). In many cases, the temperature is tightly constrained, so that in fact there are only two dimensional 
parameters that can be chosen at will. Hence, if an isothermal system satisfies a linewidth-size relation and is either 
self-gravitating or subject to ambipolar diffusion, then its velocity scale Cg oc T^/^ is set by the assumed temperature 
and its size and mean density are determined by dimensionless parameters; in this case a given simulation applies to 
only a single set of parameters describing the box. 

We return to the simplest case in which there is neither self-gravity nor ambipolar diffusion. Interstellar densities 
are often given in terms of number densities; we use the the density of hydrogen nuclei, nn = p/ fJ-H, where /in is the 
mass per hydrogen nucleus (= 2.34 x 10^^^ g for cosmic abundances). Numerically, we have for the mass, fiow time, 
and column density of the box, 






--P^l 



34.6nH,34,p. 



Mq, 



Mcs 



= 5.19x 10^ 



^0. 



pc 



.MTi 



1/2 



yr, 



(A2) 
(A3) 
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iVH = nH4 = 3.09x 102inH,34,pc cni-2, (A4) 

where Ti = T/{10 K) and c^ — 0.188 Tj^ km s^^ for molecular gas with cosmic abundances. The column density 
corresponds to a surface density 

E = 2.34 X 10"^A^H,2i g cm-2 = 11.2iVH,2i M© pc-^ (A5) 

where A^h,21 = A^h/(10^^ H cm^^). The visual extinction corresponding to this column is Ay = -/Vh.21'5 mag, 
where S is the ratio of the ext inction per unit mass to the Galactic value. The magnetic field is given by 
(jOstriker. Gammie. fc Stondll999t note that their (3 is half the normal value) 



A.n.^(^)'''=4.56(^)''' mG, (A6) 



= 3.23 (riH, 3^1)^/2-^ fiG. (A7) 

■Ma 

Scaling Relations for MHD Simulations of Turbulent Boxes with Self Gravity 

As discussed above, self-gravity introduces an additional dimensionless parameter into a simulation and therefore 
reduces the number of independent dimensional parameters by one. For the case in which the ionization scales as 

1/9 

Xi oc fijj , this reduction is the same as that due to the inclusion of ambipolar diffusion fi j2.4p : that is, for such an 
ionization law, simulations with both self gravity and ambipolar diffusion obey the same scaling relations as simulations 
with only one of these processes. In this section, we first summarize the dimensionless parameters used to characterize 
turbulent simulations with self gravity. We then describe variants of the Jeans mass that take into account turbulent 
motions. Finally, scaling laws for self-gravitating systems are given for turbulent boxes {ca = cy = 1 — see egs 1211 and 
[22|) . With two dimensional parameters specified — the strength of self gravity and the temperature — there is still one 
free dimensional parameter; as a result, self-gravitating, magnetized turbulent boxes are scale free. 

Dimensionless parameters 

There are several equivalent dimensionless parameters that can describe the effects of self-gravity. One is the ratio 
of the mass to the characteristic mass of a self-gravitating cloud, c^^/{G^pY^^: 

This parameter is related to the mass, length and sound speed by 

In terms of the free-fall velocity, Wff = {GMq/2RY/'^, we have /io = (wff/cs)'^. 
Another parameter describing the effects of self-gravity is the number of Jeans lengths in the box (e.g. 

lOstriker. Gammie. fc Ston"3ll999[ ). The typical Jeans length in the box is Aj = (7rc^/Gp)^/^, so the number of Jeans 
lengths in the box is 

I) 1/3 

n.T - 1^ = ^. (AlO) 

Aj V"" 

The Jeans mass for the box is usually defined as Mj = pAj, so that the number of Jeans masses in the box is 
Afo/Mj = (4/Aj)3 = n]. The corresponding value of iii& ^l3= M3{G^pfl'^ /cl = tt^/^ ~ 5.57. 

A third parameter describing the effects of self-gravity is the virial parameter, which for a spherical cloud of radius 
i?is 

"- = GM, (^^^) 

(jBertoldi fc McKeill992f) : here a^ = cl +cr^t is the total ID velocity dispersion. Self gravity is important for avir — 1 
and is unimportant for avir ^ 1- By contrast, fiQ and nj can have arbitrary values > 1 when self gravity is important. 
In further contrast to fiQ and nj, the effects of bulk kinetic energy as well as thermal energy are included in avir- For 
gas in a box, we define avir by replacing R by io/2: 

which is the same as equation (J19p in the text. There is a complication here, since Ofvir is defined with respect to the 
total velocity dispersion, a, whereas the linewidth-size relation depends only on the non-thermal velocity dispersion, 
(Tnt- For j\4^ 3> 1, there is no problem, since the two velocity dispersions are nearly the same. Relations involving the 
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Mach number that do not depend on the hnewidth-size relation can be extended to low Mach numbers by redefining 
Ai as [3(1 + <7^j/c^)]^/^; otherwise, such relations are restricted to A^^ S> 1. Bearing this in mind, we note that 
equation (|A9[) implies that avir is related to the other two parameters by 



6 \^l^y 67r V n^, 

iVazguez-Semadeni et alj (|2008[ ) derived a similar expression for spherical clouds (their a = |avir and their nj is R/Xj, 
which is half the value we use ). Th e virial parameter is also related to the ratio of the Jeans length to the sonic length 
(! j4.2.1l) : Equations ([55]) and (jA12p give the relation between the virial parameter and the sonic length as 

^ "^^ (A14) 



™ 6 Gp7W2£2 
for g = i . It follows that the ratio of the Jeans length to the sonic length is 

The parameters that describe the effects of self-gravity determine the ratio of the flow time to the free-fall time, 
which is 

/ o \ 1/2 

^ff = ^ ^ 1-37 X IO'^h's' yr- (A16) 



Relative to the free-fall time, the flow time is 

t±^lM4^^3.27^ = ^. (A17) 

ts M M ^1/2 ^ ' 

vir 

Variants of the Jeans length and Jeans mass 

We can define both large-scale and small-scale variants of the Jeans length and Jeans mass. On large scales, the 
density is close to the mean, p, but the velocity dispersion is cr = {a^^ + clY^'^ = A^Cs/-\/3. We therefore define the 
"turbulent" variants of the Jeans length and Jeans mass by replacing the sound speed Cs with the velocity dispersion 

^J.turb = pAj,turb= (^) Afj. (A19) 



The corresponding dimensionless quantities are 



Mo _ /V3 



3 



Aj,turb \M 

Both /iturb. and nj^turb are of order unity when the virial parameter is: 

5 5 



A'turb, 0=^3/^^3/2 -1/2)- [m) ^0' (^20) 

nj,turb^T^^=('^Vj- (A21) 



(A22) 



2Mtl.o 27rn^^,,,,- 



On small scales, however, the velocity dispersion is about equal to the sound speed, c^, whereas the de nsity can vary 
over orders of magnitude. In star- forming cores, the typical pressure is the mean turbulent pressure paf^^ ()Padoanlll995l : 
IKrumholz fc McKed 120051 ): for an isothermal gas, this corresponds to a density Pcovc = {-M.^ /3)p. We now introduce 
another variant of the Jeans length, the "core Jeans length," Aj^core, h^ which the velocity dispersion and density are 
those expected in star-forming cores, 



.2 \ 1/2 / _„2 \ 1/2 



For supersonic flows, the core Jeans length is indeed small compared to the turbulent Jeans length, Aj_ core = 
(3/A^2)Aj^turb; since it measures the effect of high pressures on thermally supported gas, whereas Aj^turb measures 



MHD Turbulence Simulations with Ambipolar Diffusion 19 

the effect of the turbulence on gas at the average density. If self-gravity is important in the turbulent box (ckvir ~ 1), 
the core Jeans length is somewhat greater than the sonic length, 

^^{^^y'=3MalC. (A24) 

The "core Jeans mass" is smaller than the normal Jeans mass and much smaller than the turbulent Jeans mass, 



Mj,eorc = PcorcAl,o„ = ^M' ^ , core = [j^j M.J. (A25) 

We expect Afj core (or perhaps the somewhat sm aller core Bonnor-Ebert mass) to be the typical mass of gravitationally 
bound cores in a turbulent cloud (|Padoanlll995[ ). 

Scaling in Terms of the Mean Density 

As discussed at the outset, the introduction of an additional physical process, such as self-gravity, reduces the number 
of independent dimensional parameters to two, which we take to be the density nu and the temperature T . Since the 
temperature has little variation in molecular clouds, there is effectively only one independent dimensional parameter, 
riH- In terms of tt-h, T and the dimensionless parameters describing the self gravity, the size and mass of the box are 
then given by 

/ y \ 1/2 / T \ ^^^ Ai f T \^^^ 

4.pc = 0.488mJ^M^^ =0.865nj^^ =0.445^7^^^ , (A26) 

jt = ^-"1 -^ = 22.3 ^^ = 3.05 ^^^ , (A27) 

ATh, 21 = 1.50Aiy'(nH, 3^1)1/2 = 2.67nj(nH,3Ti)i/^ = 1.37 ("l^Si^") M. (A28) 

Note that the column density is directly proportional to the square root of the thermal pressure in the first two cases, 
and to the square root of the turbulent pressure, pv^^^ ex M^uyiT, in the third case; this is to be expected, since the 
pressure in a self-gravitating system is proportional to GE^. The flow time is given by equations (|A16p and (|A17p . 
the magnetic field by equation (jA7|) and the ratio of the mass to the magnetic cri tical mass, M$.ni by equatiq i i (I27II . 
The scaling relation between t he size and density in terms of nj has been given bv lOstriker. Gammie. fc Stond ()1999f ) 
and IVazguez-Semadeni et all (|2005aD : for the mass in terms of the size and nj bv lTillev fc Pudritd (12004 ) although 



they have a different numerical coefficient than implied by the above relations; and bv iKlessen et al.l (|2000l ) for both 
the mass and the size in terms of the density for the particular case they consider, which has nj = 4. 
The Jeans mass and the Bonnor-Ebert mass are 



Mj = 22.3 4- Mq, Mbe = 4.74 {^] Mq, (A29) 





where for the Jeans mass we have assumed that nn is the mean density in the ambient medium, and for the Bonnor- 
Ebert mass we have assumed that nn is the density at the surface of the Bonnor-Ebert sphere. Note that Mj and 
Mbe can be expressed in terms of the surface density of the box as, for example, by 

Mj = 3.36 ^^^-^ Mq, Mbe = 0.713 ^^ ^ Mq. (A30) 

\ JVh22 I \ jVh,22 I 

We have changed the normalization of the column density so as to yield values of the Jeans mass and Bonnor-Ebert 
mass comparable to observed values for AfH,22 = Nyi/{1Q cm~^) ~ 1. The core values are smaller by a factor yji/ M. 
In this case, it is convenient to express the results in terms of the virial parameter of the box, 

Mj,corc = 5.31 ( -—^ ) Mq, Mbe, core = 1.13 ( -^ ) Mq. (A31) 

V "vir ^H, 22 / V "vir ^H, 22 / 

Note that in both cases, the critical masses have no explicit dependence on the Mach number. 

The basic conclusion is that MHD simulations of of self-gravitating, turbulent boxes are scale free; even with the 
temperature fixed, there is one free parameter, such as the box size or the density, that can be chosen arbitrarily. 
Similarly, MHD simulations with ambipolar diffusion are scale free, as we have seen in fJH If the ionization scales as 

l/rijj , then -Rad and ctvir are directly related fi j2.4p and simulations with both ambipolar diffusion and self gravity 
are also scale free. 
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Scaling with the Linewidth-Size Relation 
Linewidth-Size Relations for Molecular Clouds and Turbulent Boxes 

Molecular gas in the Galaxy exhibits a linewidth-size relation in which the velocity dispersion of the gas increases 
as a power of the physical dimension of the region (Larson 1981). Heyer & Brunt (2004) have shown that this applies 
within GMCs as well as among different GMCs. The data are consistent with the relation for the ID velocity dispersion 

fTnt = <JpcRlc, (A32) 

with cTpc = 0.72 km s"^ and 9 = 5 (jSolomon et al.l 119871: iMcKee fc Ostrikeiil2007[) : for these parameters, the Mach 

number is A^ = 6.63(i?pc/Ti)-'^/^. Such a relation appears to be satisfied by most molecular gas in the Galaxy; 
for example, Heyer fc Brunt (2004) find q = 0.49 ± 0.15 within individual molecular clouds in their sample, and 
iFalgarone et al.l (|2009f) find that the combined data from several different surveys shows a clear linewidth-size relation, 
although they do not give a fit. The parameters have different values in regions of high-mass star formation, however: 
in such regions, cTpc ~ a few km s~^, different regions of high- mass star formation do not have line widths that increase 
as i?^/^, and it is not known how the velocity dispersion within individual regions scales with size (Plume et al 1997). 
Nonethele ss, equatio n (IA32I) ap pears to be satisfi e d in t hese regions to within an order of magnitude. 

Mousch ovias fc P saltis (1995[) and iHever et al.l (l2009f) hav e proposed variants of the linewidth-size relation for gas 
that is gravitationally bound. Fal garone fc McKed ([2010) have reconciled these virialized linewidth-size relations 
with the classical turbulence- dominated linewidth-size relation. Recall that the surface density of the cloud is So = 
Mq / (cAio) , where ca = (1, 7''/4) for a box and a spherical cloud, respectively. The virial parameter (eq. IA12[) is then 

where we have used the identity symbol to emphasize that this follows directly from the definitions of the quantities 
involved; there is no physics in this relation. Solving this relation for the velocity dispersion gives 



5 (^ ' '^"^"■^o-^o 



1/2 

(A34) 



where i?o = ^o/2 is the cloud radius. The factor in parentheses with ca is unity for a spherical cloud. This relation 
is superficially like the linewidth-size relation in equation (IA32[) . but it is quite different: First, the exponent in the 
linewidth-size relation, g ~ ^ 1 follows from observation, whereas that in equation (jA34[) is ^ by definition, and second 
the coefficient in relation (jA34| depends on the column density. iHever et al.l ()2009f ) inserted the physics into this 
relation by noting that gravitationally bound clouds have avir — 1- For spherical clouds, they then found 

"TT 



—GT,oRo 
5 



1/2 

(A35) 



= 0.55(7VH.22i?pc)'/' kms-i (A36) 

which we t erm the virialized line width-size relation. They obtained a sample of bound clouds by combining the ^^CO 
data on the lSolomon et all (jl987[) molecular clouds with data from the higher resolution ^'^CO data on these clouds from 
the Galactic Ring Survey (REF). Over a range of surface densities 10 Mq pc~^ £ S < 10'^ M p, pc~^, equation (IA36I) 
describes the data well, after a somewhat uncertain correction is made for the cloud masses. iMouschovias fc PsaltisI 
(119951) previously foun d an a nalogous relation for magnetized clouds with i?^$,o hi place of N^ (see eq. [?7| . 

IFalgarone fc McKed (J201C1D concluded that non-self gravitating interstellar gas obeys the turbulence-dominated 
linewidth-size relation given by equation (|A32[) . whereas self-gravitating gas satisfies the virialized linewidth-size rela- 
tion given by equation (IA35p . There is a critical surface density that defines the boundary between the turbulent and 
virialized cases: Equating the velocity dispersions for the two cases gives 

= 192(^^)a;/ Mopc^, (A38) 

where cr*^ = crpc/(0.72 km s^^) is normalized to the standard Galactic value (see eq. [S5]) and ca = vr/4 for spherical 
clouds. The corresponding column density is 

^Lws = 1.71 X 10^2 {^\ a;/ cm-2. (A39) 

It should be noted that these value s for Sl ws and A'lws for spherical clouds are comparable to the mean values 
for Galactic GMCs found bv .Solomon et al.l ^1987.) . and about twice the mean values found by IHever et al.l pOQ9.) . 
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Regions with E <C Slws are dominated by interstellar turbulence, whereas those with E ^ Elws are dominated by 
self-gravity and are decoupled from the turbulent cascade in the interstellar medium. The high surface densities of 
regions of high-mas s star formation thus naturally lead to the high velocity dispersions observed there by, for example, 
iPlume et alJ (|1997l ) . The maximum density for a cloud satisfying the turlDulent linewidth-size relation is 

where we used equation ([M)) to eliminate t and the relations M (x caN^I^ — cynH^^ to cover the different geometries. 
The maximum mass for a cloud satisfying the turbulent linewidth-size relation is 

Mums = 0.311 I — ^ J M©, (A41) 

which is independent of geometry. 

Turbulent Boxes with the Turbulence- Dominated Linewidth-Size Relation {Nh < Nuws, ca = cy = 1) 

Here we determine the scaling relations for turbulent boxes that satisfy the turbulence-dominated linewidth-size 
relation 

cr„t = 0.72o-;^i?J/2 kms-i. (A42) 

The Mach number A^ and box size Iq are related by equation ([5^ . As a result, the properties of the turbulent box 
are given by: 



:2.36xl05(— ^ yr, (A44) 



Sc 



Mo = 3.25xlO-^( -^'_yn Mo, (A45) 






iVH = 1.40xlO^" pj^f^^^M cm-2. (A46) 

For simulations that include self- gravity but have surface densities less than the critical one (E < Elws), the column 
density is most simply expressed in terms of the linewidth-size parameter and the virial parameter using equations 
(|A33l) and (lA43| . 

N^ = U^\-^ = l.MxlQ^^{''-^\ cm-2, (A47) 

4 ylpcy G/^navir y avir y 

which corresponds to iVn = iVLWs/ctvir for ca = 1 in equation (IA39p . The scaling for the density can be expres sed in 
terms of the linewidth-size parameter a* and a parameter describing the self gravity with the aid of equations (jA43p 
and (|X26)) , 

(2/3*4\ /2*4\ / *4\ 

f^!L^ ^36lK^U96 ^P- . (A48) 

Comparison with equation (jA40p for a box geometry (cy = 1) shows that njj = "Lws/c^virj where avir ^ 1- Similarly, 
equation (jA27|l implies 

so that Mq — MuMs/oivir, where Mlws is given in equation (jA4ip . The turbulence-dominated linewidth-size relation 
does not apply for densities exceeding tilwSj correspondi ng to c olumn densities N^i > A^lws and masses M > AfLvvSi 
and for that case the scaling is given by the relations in i 3A.2.3l below. 

When expressed in terms of the linewidth-size relation, the core values for the Jeans mass and Bonnor-Ebert mass 
(eq. IA25P are independent of the Mach number, 

Mj,core = 3.96^^i^ M0, MBE,core = 0.84^^i^ Mq. (A50) 
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The core Bonnor-Ebert mass is comparable to the typical mass of observed stars, pa rticularly if allowance is m ade for 
the fact that only a fraction of the core mass is incorporated into the final star (e.g. iMatzner fc McKeil2000f) . These 
relations can be expressed in terms of the sonic length instead of a* by using equation ((39)) . 

In applying these relations to simulations with drive n tur bulence, it must be kept in mind that the driving generally 
results in deviations from the linewidth-size relation (|A32p on the driving scale. Using the linewidth-size relation to 
relate simulations to actual systems is therefore best done for cases in which the driving is restricted to large scales; 
the simulations discussed in the text satisfy this constraint since they are driven over a narrow range of wavenumbers 
at the largest scale, 1 < fe < kd, with the driving wavenumber kd = 2. (Here fc is a dimensionless wavenumber that 
is related to the physical wavenumber fcphys by fc = fcphys^o/27r; the minimum possible wavenumber is fc = 1 and the 
maximum is Afg/2, where Afg is the number of grid cells in each side of the box.) When the turbulence is driven, one 
must distinguish between the linewidth-size relation applied to the entire box, and the linewidth-size relation inside 
the box. In our simulations, the mean Mach number is approximately constant over the range 1 < fc < fc^ = 2. We 
have chosen to use the full size of the box in relating our simulations to clouds: io — 2R, where R is the radius of the 
cloud or of a region inside the cloud. However, in determining properties inside the cloud, such as the sonic length, it 
is necessary to allow for the fact that the internal linewidth-size relation is normalized approximately to the driving 
scale. As a result, for a oc i^^"^, the sonic length is is — {^o / kd) / M"^ rather than Iq/M'^ (see eq. | 



General Scaling with the Lmewidth-Size Relation 

The virialized linewidth-size relation fo llows f rom assuming that the virial parameter is unity, so the scaling relations 
for this case are given by the results in JjA. 1.31 with avir = 1. To cover both the turbulence-dominated and virialized 
cases, note that equation (JA46I) shows that A^h oc tt-h for the turbulence-dominated case (avir ^1), corresponding to 

"■H ^ ^LWS: whereas equation (jA28[) shows that A'h oc n^ for avir — 1, corresponding to fin ^ nims- As a result, we 
have 

/ - \ 1/2" 



A^H = A^LWS min 



"LWS V "LWS 



(A51) 



for the turbulence-dominated and virialized cases, respectively, as can be verified by direct substitution using equations 
(PS8)) . (|X39|l and (1X40)) . Similarly one can show that 



Mo = Mlws min 



/ - N 1/2' 



"LWS V "-H 



(A52) 



with the aid of equations ()A27I) , ()A40I) and ()A41I) . Note that Mlws is the maximum possible mass for a cloud with 
a given velocity dispersion, a ex A^T^/^, and linewidth-size coefficient, a* (jFalgarone fc M cKcc 2O1O0 . Furthermore, 
for a simulation with a given Mach number, Mlws decreases as a* increases (eq. IA41[) : as a result, Mlws is also the 
maximum mass of a cloud with a linewidth above the linewidth-size relation. The size of the simulation box is 

1/21 



0.0455 ^ 



a* ^ 



1, p^ 



pc (A53) 



based on equation ()A26|) . It must be borne in mind that these equations are based on the mean linewidth-size relation; 
for a given size and/or surface density, the velocity dispersion can vary by a factor of a few. Thus the virialized 
linewidth-size relation, which applies to regions with high column densities by definition, correspondingly applies to 
regions of high density but with sizes and masses that decrease as the column density increases. 

Code Units 

Numerical codes are generally written in dimensionless form, with masses, lengths and times written in terms of 
code units, M = M/Mcode, ^ = ^/4ode, and i = t/tcode- The code units can be adjusted to fit the problem being 
simulated. The properties of the box in code units, Mg = -^o/-^code and ^o — ^o/^codc: can be selected arbitrarily 
prior to the simulation (e.g., Mq = 8 and i'o = 2), as can the normalized sound speed, c^ — Cgtcodc/ ^mdc- If there 
are Mg grid cells in each side of the simulation box, the grid size is lS.t — £o/Afg. For stationary gas, the time step is 
At = CA£/cs, where C is the Courant number. 

The code unit for length is given by ^codc = ^o /^o, where £o is given by equation (|A43[) if the typical Galactic 
linewidth-size relation is adopted, by equation ()A26I) for a self-gravitating gas, and by equation (jTT)) for a gas undergoing 
ambipolar diffusion. The corresponding code unit for time is given by tcodc = Csicodc/cs- The code unit for mass is 
given by Afcodc = Mq/Mq, where Mq is given by equation ()A27|) for a self-gravitating gas, by equation (|A49I) for 
a self-gravitating gas that obeys the linewidth-size relation, and by equation (1491) for a gas undergoing ambipolar 
diffusion. 

The gravitational constant in the code is 

tcodctld. ^ (GMo\ hB- = ^0^'/°^' 



^codc V 4c2 ; Mo Mo 
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from equation (jA9p . Including the Heavy Ion Approximation (xi 
constant in the code is 



TZXi and 7ad oc jab/T^-), the ambipolar diffusion 



7AD 



lAuMcodct 



codc'^codc 



TZi 



3 
code 



2i?AD(4) 

nx^M|3 



Mn' 



(A55) 



where the second step follows from equation ([8|). [Keep in mind that £o, Cs, and Mq are arbitrary; in the simulations 

described in the text, we have taken £q = 2, Cg = 0.1 and p — Mq/£q = 1, so that 7ad — 0-^RAD{£o)/T^Xi-^P-] So 
long as the Heavy Ion Approximation is valid, the outcome of a simulation is independent of the value of TZ since 7ad 
always enters in combination with Xi = T^Xi- 

We can now address the issue of scaling in AD simulations when the physical ionization is specified by, for example, 
the value of x*- In carrying out a simulation, the ionization in code units, Xi — ^Xii must be specified; hence, the 
physical ionization Xi oc Xi '^ TZ~^ . As noted above, the results of a simulation are independent of 7?. so long as the 
Heavy Ion Approximation is valid. Thus, a single simulation provides the results for a family of problems with different 
degrees of ionization but the same values of i?AD(^o) and c^; as a result, we can use a single simulation to treat the 
physically plausible range of ionizations for a given value of i?AD(^o): a-s discussed in §5.11 
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Fig. 1. — Clumps' AD Reynolds number, -Rad.o versus ion Alfvcn Mach number squared, A^Ai^i for models m3e2r-l(eircles), 
m3c2rl (squares), and m3e2r3 (diamonds) at the end of the simulation. The straight line shows -Rad.c = AIai and the solid symbols 
indicate the values of Rad and AiAi for the whole box. 
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Fig. 2. — Clump mass functions for AD models of i?AD = 1200 with grid sizes of 256^ (dashed line) and 512^ (m3c2r3, solid line). The 
sonic wave number, ka, and the minimum scale of the inertial range, fcin.maxi for model m3c2r3 are plotted as vertical solid lines; fcin.max 
for the 256'^ model is plotted as a vertical dashed line. The clump wavenumber kc based on equation 1 140 II is plotted at the top of the figure 
for reference; note that k increases to the left. See 54.2.1 for discussion. 
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Fig. 3. — 3D spatial distribution of clumps, identified by CLUMPFIND with minimum mean radius of 6 cells, from model m3c2rl. 
Different gray-scale shadings (different colors in the online version of the paper) are used only to visually separate overlapping individual 
clumps. 
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Fig. 4.— Clump mass functions for models m3c2r-l [Rad(^o) = 0.12; solid line], m3c2rl [HadC^o) = 12], and m3i (ideal MHD). The 
dashed line and dot-dashed line show the best fitting higher-mass slope Fflt for m odels m3c2rl and m3i, respectively. Model m3c2r-l, which 
has the strongest ambipolar diffusion, has the steepest higher-mass slope (see ^4.2.21 for discussion); Model m3i is very similar to model 
m3c2r3 [-Rad(^o) = 1200], which is not shown. All clumps with radius larger than 3 cells are included in the plot. 
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Fig . 5 . — Convergence study for the mass-to-flux ratio of the clumps . Values of normalized mass-to- flux ratios of clumps in the -Rad = 1 200 
model are plotted against normalized clump mass, for resolutions of 256 (blue cirlces) and 512 (red squares; model m3c2r3). By plotting 
(M*,c/M#,o)/(A^c/Afo)^'^ versus Mc/Mq, data points are projected horizontally for easy visual comparison. The straight line shows the 
mean of the normalized mass-to-flux ratio for elliptical clumps using the mean density and B-field of the clumps from model m3c2r3. See 
54.3.1 for discussion. 
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Fig. 6. — Normalized mass-to-flux ratios, (M*,c/M4>,o)/(A^c/Afo)"'^' '*, of clumps in models m3c2r-l (blue circles), m3c2rl (black crosses), 
and m3c2r3 (red triangles) plotted versus normalized clump mass Mc/Mq. The mean values of {fJ.^,c/fJ-^,o/{Mc/Mo)^'^) for the three 
models are plotted as the horizontal lines (-Rad(^o) = 0.12 blue solid, -Rad(^o) = 12 black dashed, and R\y>{£o) = 1200 red dot-dashed). 
The model with the largest value of -Rad has a slightly lower average mass-to-flux ratio and a smaller dispersion of the mass-to-flux ratios. 
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Fig. 7. — Other normalized physical properties of clumps as functions of the AD Reynolds number of the clumps, -Rad, o for the models 
m3c2r-l (on the left) and m3c2r3 (on the right): (a) radius, rc/^Oi (b) ion density; {pi)^/p\ (c) neutral density, {pn)i,/p', (d) magnetic 
energy density, Ub,c/Ub,o\ (g) clump mass, Mc/Mo; and (f) ionization mass fraction, Xi.c/Xi.O- See i|4.4l for discussion. 
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TABLE 1 

AD Reynolds Number Rad for Observed Molecular Clumps (|CrutcherII1999I ) 



■CEHH 3 IS?TIi 7? M Ma 31 — ^^ 

(Ha cm-^') (pc) (K) 



W3 0H 


0.07 


6.8 


0.02 


1.9 


0.3 


100 


3.0 


DR 21 OHl 


0.21 


6.3 


0.05 


4 


1.3 


50 


37.3 


SgrB2 


0.0008 


3.4 


22 


22 


0.4 


70 


10.3 


M17SW 


0.008 


4.5 


1 


7 


0.5 


50 


6.3 


W3 (main) 


0.13 


5.5 


0.12 


4.8 


1.2 


60 


24.1 


S106 


0.04 


5.3 


0.07 


3.6 


0.5 


30 


3.7 


DR 21 OH2 


0.41 


6 


0.05 


4 


1.8 


50 


51.5 


OMC-1 


0.65 


5.9 


0.05 


1.7 


1 


100 


21.9 


NGC 2024 


0.35 


5 


0.2 


3.7 


1.6 


25 


72.7 


S88B 


0.056 


3.8 


0.7 


5.9 


1 


40 


12.9 


Bl 


0.17 


4 


0.2 


3.6 


1.1 


12 


15.7 


W49 B 


0.024 


3 


1 


5.9 


0.6 


10 


6.3 


W22 


0.033 


3 


4 


3.5 


0.5 


10 


20.5 


W40 


0.027 


2.7 


5 


10 


1.2 


10 


42.4 


p Oph 1 


0.42 


3.2 


0.8 


3.5 


1.6 


25 


41.6 


OMCN-4 


>0.47 


6 


0.03 


2.9 


>1.4 


35 


>30.7 


Tau G 


>0.042 


3 


1 


5.1 


>0.7 


10 


>9.5 


L183 


>0.052 


3.1 


0.3 


2.4 


>0.4 


10 


>1.9 


L1647 


>0.047 


3 


3 


9 


>1.4 


10 


>56.4 


pOph2 


>0.14 


3 


0.9 


3.2 


>0.8 


25 


>11.3 


TMC-1 


>0.063 


3 


1.9 


5.9 


>1 


10 


>31.4 


L1495 W 


>0.063 


3 


0.9 


3.9 


>0.7 


10 


>9.8 


L134 


>0.14 


3.2 


0.3 


2.7 


>0.7 


10 


>6.3 


TMC-IC 


>1.3 


4 


0.2 


2 


>1.6 


10 


>73.0 


L1521 


>0.13 


3 


1.2 


3.9 


>1 


10 


>27.0 


L889 


>0.28 


3 


2.4 


7.3 


>2.7 


13 


>191.1 


Tau 16 


>0.22 


3 


1.2 


3.9 


>1.3 


10 


>45.7 



" i?AD computed using equation II15I I 
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TABLE 2 

Model Parameters and Regimes of AD 



Model" 


7AD 


Rad(4) 


RADim 


Regime of AD 


m3c2r-l 


4 


0.12 


0.076 


III 


m3c2r0 


40 


1.2 


0.70 


II ~ III 


m3c2rl 


400 


12 


10.1 


II 


m3c2r2 


4000 


120 


103.2 


II 


m3c2r3 


40000 


1200 


1022 


I 


m3i 


oo 


oo 




I 



" Models are labeled as "rnxcyrn," where x is the thermal Mach number. 
MHD. Model m3c2r0 is the same as model m3c2h in LMKF. 
' /?AD from models using time-dependent ionization (see §3). 
"^ Root mean squared (rms) values. 



y = I logXioli S'lid n = log(R\Y^(£)/1.2). Model "m3i" is an ideal 
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TABLE 3 

Comparison of Clump Properties in Models with Different Rao 



Model 


m3c2r-l 


m3c2rl 


m3c2r3 


RAD(io) 

n„„(fc)" 
(m*,c>/m|,o 

(Rc./Rc,±) 

(rcXcell)- 

(Nc) (dc > 12 cells)-'' 


0.12 
1.96 ±0.02 

0.122 ±0.004 

0.056 

6.85 ±0.24 

1.13 ±0.02 

9.4 ±1.1 

698 


12 

1.89 ±0.03 

-1.43 ±0.10 

0.120 ±0.005 

0.074 

5.42 ±0.23 

1.15 ±0.02 

12 ±1.5 

434 


1200 

1.48 ± 0.05 

-1.22 ±0.11 

0.111 ±0.003 

0.039 

3.19 ±0.09 

2.00 ± 0.05 

13.5 ±1.7 

349 



° Velocity power spectral index of neutral component. 

' Slope of the CIMF in the inertial range. The data for model m3c2r-l do not have a single power law over the inertial range. 

"^ Clump mass-to-flux ratio normalized by that of the whole box. 

'* Dispersion of clump mass-to-fiux ratio normalized by that of the whole box. 

* Mean radius of clumps in units of number of cells. 

f Mean number of clumps with diameter larger than 12 cells. 



